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I.  Summary 

We  consider  a  secondary  multi-hop  cognitive  radio  network  with  decentralized 
control  that  operates  cognitively  to  coexist  with  primary  users.  We  propose  a 
new  spread-spectrum  management  paradigm,  in  which,  unlike  mainstream  dy¬ 
namic  spectrum  access  research,  digital  waveforms  are  designed  to  occupy  the 
entire  available  spectrum,  and  to  adaptively  track  the  interference  profile  at 
the  receiver  to  maximize  the  link  capacity  while  avoiding  interference  to  pri¬ 
mary  users.  In  this  context,  we  study  the  problem  of  maximizing  the  network 
throughput  of  a  multi-hop  network  through  joint  routing  and  spread-spectrum 
channelization.  We  first  propose  a  centralized  formulation  of  the  network  con¬ 
trol  problem.  We  then  propose  an  algorithm  that  can  be  seen  as  a  distributed 
localized  approximation  of  the  throughput-maximizing  policy.  We  refer  to  the 
proposed  jointly-designed  routing  and  code-division  channelization  algorithm 
as  ROCH  (Routing  and  cOde-division  CHannelization).  Specifically,  power  and 
spreading  code  are  jointly  selected  to  maximize  the  pre-detection  secondary 
SINR  while  providing  quality  of  service  guarantees  to  on-going  primary  and 
secondary  transmissions,  while  the  routing  algorithm  dynamically  selects  re¬ 
lays  based  on  the  network  traffic  dynamics  and  on  the  achievable  data  rates 
on  different  secondary  links.  We  study  the  throughput  and  delay  performance 
of  ROCH  through  a  extensive  simulation  experiments,  which  demonstrate  the 
appeal  of  the  proposed  framework  through  significant  performance  gains  com¬ 
pared  to  baseline  solutions. 

Existing  commercial  wireless  systems  are  mostly  hardware-based,  and  rely 
on  closed  and  mostly  inflexible  designs  and  architectures.  Moreover,  while  in 
recent  years  there  have  been  significant  algorithmic  developments  in  cross-layer 
network  adaptation  and  resource  allocation,  in  most  cases  existing  architectures 
are  not  flexible  enough  to  incorporate  these  advancements  other  than  through  ad 
hoc  tweaks.  To  address  these  limitations,  we  propose  RcUBe  (Real-time  Recon- 
figurablc  Radio),  a  novel  radio  framework  that  offers  real-time  reconfigurability 
and  self-optimization  capabilities  at  the  PHY,  MAC,  and  network  layers  of  the 
protocol  stack.  Unlike  state-of-the-art  solutions  designed  to  provide  flexibility 
at  a  specific  layer  of  the  protocol  stack,  RcUBe  is  a  flexible  architectural  ab¬ 
straction  spanning  the  three  lowest  layers  of  the  protocol  stack.  RcUBe  provides 
abstractions  to  define  in  a  structured  way  sophisticated  decision-making  algo¬ 
rithms  and  applying  the  corresponding  reconfigurations  in  real-time.  RcUBe 
provides  these  features  through  a  design  structured  into  four  planes,  namely 
decision ,  control ,  data ,  and  register  plane.  The  decision  plane  consists  of  user- 
defined  decision  algorithms  that  provide  the  control  logic  for  real-time  protocol 
optimization,  adaptation,  and  reconfiguration.  In  the  control  plane,  RcUBe 
stores  the  required  logic  for  routing  and  MAC  protocol  execution  as  well  as 
data  plane  management.  Data  processing  takes  place  in  the  data  plane,  while 
the  register  plane  stores  and  manages  access  to  system  parameters  and  environ¬ 
mental  information.  We  illustrate  the  unique  features  and  the  flexible  design 
of  RcUBe  by  implementing  a  cross-layer  architecture  that  performs  on-the-fly 
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reconfigurations  based  upon  decisions  taken  in  real  time  by  a  decentralized  op¬ 
timization  algorithm. 

Finally,  we  design  and  implement  a  wireless  software  defined  radio  (SDR) 
testbed  for  cognitive  channelization  in  multipath  fading  environments.  Pri¬ 
mary  and  secondary  users  coexist  and  utilize  the  same  spectrum  concurrently. 
Cognitive  channelization  is  achieved  by  adaptive  waveform  assignment  to  sec¬ 
ondary  users  that  assumes  unknown  primary  users  and  maximizes  the  signal- 
to-interference-plus-noise  ratio  (SINR)  at  the  output  of  the  secondary  receiver, 
while  at  the  same  time  maintains  interference  to  primary  users  below  a  prespec¬ 
ified  threshold.  We  deployed  four  commercial,  low-cost  SDR  transceivers  in  an 
indoors  lab  environment.  To  address  challenges  related  to  low-cost,  off-the-shelf 
testbed  components  and  in  particular  to  unstable  clock  oscillators  that  intro¬ 
duce  carrier  frequency  offsets  (CFOs)  between  any  transmitter-receiver  pair, 
we  propose  a  subspace  based  technique  for  joint  blind  estimation  of  the  multi- 
path  channel  coefficients  and  the  introduced  CFOs.  The  theoretical  concepts 
of  adaptive  cognitive  channelization  in  the  presence  of  unknown  primary  users 
and  multipath  fading  are  experimentally  demonstrated  and  validated  in  terms 
of  instantaneous  SINR  and  bit-error-rate  (BER). 
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II. 


Introduction 


A.  All-spectrum  Cognitive  Networking  through  Joint  Dis¬ 
tributed  Channelization  and  Routing 

Cognitive  radio  networks  [l]-[2]  have  emerged  as  a  promising  technology  to 
improve  the  utilization  efficiency  of  the  existing  radio  spectrum.  Mainstream 
cognitive  radio  proposals  focus  on  opportunistic  access  to  the  licensed  spectrum 
where  the  primary  users  of  the  band  are  known  a  priori  and  this  knowledge 
can  be  utilized  to  detect  if  the  band  is  occupied  by  the  known  signal  pattern. 
Quite  the  opposite,  in  the  unlicensed  band  there  are  potentially  many  unco¬ 
ordinated  devices,  and  their  signal  waveforms  and  activation  statistics  are  in 
general  unknown.  Moreover,  in  cognitive  radio  networks  with  multi-hop  com¬ 
munication  requirements,  spectrum  occupancy  is  location-dependent,  and  the 
receiver  interference  profile  may,  thus,  vary  at  each  relay  node. 

We  therefore  propose  a  new  spread-spectrum  management  paradigm,  in 
which  waveforms  are  designed  to  occupy  the  entire  available  spectrum  with¬ 
out  generating  harmful  interference  to  active  primary  or  secondary  users.  In 
this  way,  the  secondary  users  share  the  licensed  spectrum  with  the  primary 
users  to  achieve  frequency  reuse.  At  the  same  time,  the  dynamic  and  location- 
dependent  nature  of  the  wireless  environment  calls  for  the  development  of  rout¬ 
ing  algorithms  that  are  aware  of  the  interference  profile  at  each  potential  relay. 

The  lack  of  established  infrastructure  and  the  wireless  channel  dynamics  im¬ 
pose  an  unprecedented  set  of  challenges  over  spread-spectrum  cognitive  ad  hoc 
networks.  First,  secondary  users  should  optimize  the  spreading  code  and  power 
to  avoid  generating  harmful  interference  to  primary  users.  The  challenges  here 
arise  from  the  assumption  that  the  spreading  codes  of  primary  users  are  un¬ 
known  to  the  secondary  users.  Second,  in  a  multi-hop  network,  the  spectrum 
environment  varies  in  time  and  space  depending  on  the  activities  of  primary 
users,  interference,  and  fading.  The  optimal  spectrum-spreading  channelization 
may  therefore  be  different  at  each  hop  in  a  multi-hop  path.  Furthermore,  as 
new  secondary  links  are  formed  and  others  vanish,  and  following  network  traf¬ 
fic  dynamics,  routing  of  data  flows  from  one  secondary  node  to  another  may 
frequently  change.  Therefore,  controlling  the  interaction  between  routing  and 
code  design  is  of  fundamental  importance. 

In  this  work,  we  explore  a  new  framework  that  captures  the  interdepen¬ 
dencies  between  spread-spectrum  channelization  and  routing.  The  through¬ 
put  optimization  is  carried  out  dynamically  by  all  secondary  transmitters  to 
continuously  adapt  to  the  changing  spectrum  environment  and  traffic  arrival 
rates.  Specifically,  a  distributed  algorithm  for  dynamic  joint  power  and  spread- 
spectrum  channelization  is  developed  to  maximize  the  pre-detection  secondary 
SINR  while  guaranteeing  the  SINR-QoS  requirements  for  on-going  transmissions 
from  primary  and  secondary  users.  The  excellent  cognitive  network  performance 
characteristics  is  demonstrated  by  simulation  studies  included  in  this  project. 

A  motivating  example.  Through  an  example,  we  will  try  to  intuitively 
highlight  the  potential  benefits  of  spread-spectrum  channelization  as  compared 
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to  traditional  frequency-division  multiplexing  (FDM)  for  dynamic  spectrum  ac¬ 
cess.  In  code-division  channelization,  secondary  users  use  spread  spectrum  sig¬ 
nals  that  occupy  the  whole  spectrum  band.  In  this  example,  16  —  bit  spreading 
codes  over  the  real  field  are  assigned  to  secondary  users.  In  FDM,  the  whole 
spectrum  band  is  divided  into  16  subbands,  and  each  primary  user  is  assigned 
with  a  specific  subband.  In  FDM  scenario,  secondary  users  can  access  a  single 
subband  only  if  the  subband  is  not  occupied  by  primary  users  or  any  other 
secondary  users.  In  the  simulation,  the  transmission  SNRs  of  primary  links  are 
all  set  equal  to  15  dB;  the  maximum  allowable  transmission  SNR  for  secondary 
links  is  set  to  15  dB.  All  the  signature  vectors  for  primary  links  are  generated 
from  a  minimum  total-squared  correlation  optimal  binary  signature  set  which 
achieves  the  Karystinos-Pados  (KP)  bound  [32]- [34]  for  each  (K,  16)  pair  of  val¬ 
ues1,  where  K  is  the  number  of  active  primary  users.  Both  of  code-division 
channelization  and  FDM  schemes  provide  the  same  minimum  SINR  guarantee 
3  dB  for  both  active  primary  users  and  secondary  users. 

Note  that  we  require  that  both  the  code-division  channelization  and  the 
FDM  schemes  provide  the  same  minimum  SINR  guarantee  for  both  active  pri¬ 
mary  users  and  secondary  users.  We  vary  the  number  of  active  primary  users 
and  evaluate  the  maximum  number  of  admissions  for  secondary  users  with  the 
two  different  schemes.  As  shown  in  Fig.  1,  it  can  be  observed  that  spread 
spectrum  management  with  code-division  channelization  allows  more  secondary 
users  to  simultaneously  access  the  spectrum  than  FDM  spectrum  management. 
For  example,  with  16  subbands  and  14  active  primary  users,  only  2  secondary 
users  can  access  the  network  with  FDM.  While  up  to  7  secondary  users  can 
access  the  network  with  carefully  designed  code-division  channelization. 

Contributions.  Within  this  context,  our  main  contributions  in  this  project 
can  be  outlined  as  follows: 

•  Uncoordinated  code-division  spectrum  management.  Unlike 
mainstream  work  on  cognitive  radio  networks,  we  consider  a  code-division 
channelization  based  secondary  network  that  operates  cognitively  to  co¬ 
exist  with  primary  users  for  infrastructure-less  cognitive  radio  ad  hoc  net¬ 
works. 

•  Distributed  joint  routing  and  code-division  channelization.  We 

formulate  a  joint  routing  and  code-division  channelization  problem.  Given 
the  centralized  nature  and  high  computational  complexity  of  the  prob¬ 
lem,  we  study  localized  algorithms  for  joint  dynamic  routing  and  spread- 
spectrum  channelization  that  are  designed  to  maximize  the  global  objec¬ 
tive  function  of  the  centralized  problem.  To  the  best  of  our  knowledge, 
no  existing  algorithm  attempts  to  control  the  routing  and  code-division 
channelization  functionalities  to  jointly  maximize  the  network  throughput. 

•  Novel  low-complexity  spread-spectrum  channelization  algorithms. 

We  study  the  problem  of  designing  a  joint  power  and  spreading  code  al- 

1When  K  <=  16,  the  KP-optimal  sequences  coincide  with  the  familiar  Walsh-Hadamard 
signature  codes. 
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Fig.  1:  Comparison  between  code-division  channelization  and  FDM. 


location  scheme  for  the  secondary  link  that  maximizes  the  pre-detection 
SINR  under  quality  of  service  constraints  for  the  secondary  users.  Us¬ 
ing  semidefinite  relaxation,  we  develop  a  novel  low-complexity  suboptimal 
solution  to  an  otherwise  NP-hard  problem  with  excellent  performance  in 
practice. 

B.  Related  Work 

Since  spectrum  occupancy  is  location-dependent,  the  spectrum  occupancy  pro¬ 
file  may  be  different  at  each  relay  node  in  a  multi-hop  end-to-end  path.  There¬ 
fore,  one  of  the  key  challenges  in  the  design  of  cognitive  ad  hoc  networks  is  to 
jointly  and  dynamically  allocate  routes  and  portions  of  the  spectrum  to  each 
node  in  the  multi-hop  network.  The  authors  in  [3]  proposed  a  routing  algorithm 
based  on  a  probabilistic  estimation  of  the  available  capacity  of  each  secondary 
link.  The  proposed  probability-based  routing  metric  relies  on  the  probability 
distribution  of  the  interference  between  primary  and  secondary  users  over  a 
given  channel.  In  [4,  5],  the  authors  proposed  a  connectivity-driven  routing 
algorithm,  where  paths  are  measured  in  terms  of  their  degree  of  connectivity 
in  a  multi-hop  cognitive  radio  network  that  is  highly  influenced  by  the  primary 
user  behavior.  Route-stability-oriented  routing  is  introduced  in  [6]  based  on 
the  concept  of  route  maintenance  cost.  The  maintenance  cost  represents  the 
effort  needed  or  penalty  paid  to  maintaining  end-to-end  connectivity.  In  [7],  a 
distributed  and  localized  algorithm  for  joint  dynamic  routing  and  spectrum  allo¬ 
cation  for  ad  hoc  cognitive  radio  networks  is  proposed.  The  proposed  algorithm 
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jointly  addresses  routing  and  spectrum  assignment  with  power  control  under  the 
physical  interference  model.  In  [8],  a  coordination  scheme  is  introduced  that  al¬ 
lows  secondary  users  to  coordinate  with  primary  users.  Specifically,  secondary 
users  offer  their  services  as  an  intermediate  relay  node  in  an  effort  to  improve 
throughput  of  primary  users  utilizing  a  802.11-based  channel  access  mechanism. 
In  return,  the  secondary  user  ‘piggy-backs’  some  of  its  own  data  while  acting 
as  a  relay.  In  [9],  the  authors  propose  the  CRP  routing  protocol  that  consid¬ 
ers  joint  spectrum-route  selection  and  service  differentiation  in  CR  routes.  The 
proposed  protocol  allows  for  two  classes  of  routes  -  class  I  routes  that  provide 
better  CR  network  performance,  while  class  II  routes  aim  to  achieve  a  higher 
measure  of  protection  for  the  PUs.  The  authors  in  [10]  proposed  a  distributed 
routing  algorithm  where  secondary  users  minimize  their  interference  to  the  pri¬ 
mary  users  while  keeping  the  delay  along  the  route  low.  The  reader  is  referred 
to  [11]  and  references  therein  for  an  excellent  survey  of  the  main  results  in  this 
area. 

Herein,  we  consider  cognitive  networks  built  around  a  primary  code-division 
multiplexing  (CDM)  system.  Unlike  traditional  frequency  division  operations 
where  cognitive  secondary  users  may  transmit  opportunistically  in  sensed  spec¬ 
trum  holes/void  only,  cognitive  code-division  users  may  operate  in  parallel  in 
frequency  and  time  to  a  primary  system  as  long  as  the  induced  spread-spectrum 
interference  remains  below  a  pre-defined  acceptable  threshold2.  Power  control 
under  an  “interference  temperature”  constraint  (total  secondary  user  distur¬ 
bance  power  over  primary  band)  was  considered  [12]  in  cognitive  code-division 
systems.  Furthermore,  a  joint  power  allocation  and  admission  control  algo¬ 
rithm  was  proposed  in  [13]  for  cognitive  CDMA  systems  to  maximize  the  en¬ 
ergy  efficiency  of  secondary  users  with  QoS  guarantee  for  primary  and  secondary 
links.  In  [14],  a  spectrum  underlay  (CDM-based)  cognitive  radio  network  (CRN) 
was  considered  and  a  low-complexity  suboptimal  algorithm,  named  interference 
constraint-aware  stepwise  maximum  interference  removal  algorithm  (I-SMIRA), 
was  proposed  to  maximize  the  number  of  admitted  secondary  users  by  optimiz¬ 
ing  power  and  transmission  rate.  However,  no  code-channel  (signature)  opti¬ 
mization  was  carried  out  for  the  secondary  users.  In  contrast,  in  [15]  a  sec¬ 
ondary  code  assignment  scheme  was  presented  to  obtain  the  binary  secondary 
signature  by  hard-limiting  the  code  sequence  that  exhibits  the  minimum  mean- 
square  cross-correlation  with  the  primary  received  signal.  The  secondary  code 
set  of  multiple  secondary  users  was  also  constructed  in  an  iterative  way.  Under 
interference-minimizing  code  assignments,  bit  rate  and  spreading  factor  adjust¬ 
ments  for  a  secondary  CDMA  system  were  considered  in  [16].  A  distributed  al¬ 
gorithm  for  resource  allocation  of  spectral  bands,  power,  and  data  rates  among 
multiple  secondary  users  for  multi-carrier  CDMA  systems  was  developed  in  [17]. 
The  problem  of  beamforming  and  power  control  for  CR  networks  aimed  at  re¬ 
stricting  the  interference  on  PUs  has  been  addressed  in  [18].  The  authors  in 
[19]  study  the  problem  of  fair  spectrum  sharing  among  all  SUs  in  underlay  CR 

2While  early  standardization  and  regulation  discussions  have  begun  [23],  no  conclusive 
“interference  temperature”  rules  and  agreements  have  been  reached  yet. 
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networks,  subject  to  certain  QoS  (in  terms  of  minimum  SINR  and  transmission 
rate)  and  outage  probability  of  constraints  with  imperfect  channel  state  infor¬ 
mation.  A  joint  admission  control  and  rate/power  allocation  for  secondary  users 
is  proposed.  In  [20],  the  authors  consider  the  scenario  where  SUs  can  adaptively 
adjust  their  transmit  directions  in  addition  to  transmit  power  according  to  avail¬ 
able  channel  information.  A  joint  power  allocation  and  phase  control  solution 
is  proposed  subject  to  interference  and  power  constraints.  In  [21],  the  authors 
study  the  optimal  scheduling  problem  with  the  objective  to  achieve  proportional 
fairness  of  the  long-term  average  transmission  rates  among  different  links  in  a 
cognitive  ad  hoc  network  with  spectrum  underlay.  Hybrid  overlay /underlay  CR 
waveforms  were  designed  to  adapt  its  spectrum  to  efficiently  exploit  both  unused 
spectrum  holes  and  underused  spectrum  bands  through  OFDM  and  MC-CDMA 
in  [22], 

The  common  theme  of  most  research  in  the  area  of  CDMA-based  CR  net¬ 
works  is  to  optimize  physical  layer  performance,  without  modeling  in  detail 
how  physical  layer  resource  allocation  interacts  with  higher  layers  of  the  pro¬ 
tocol  stack  to  improve  network  performance  metrics.  In  comparison  with  work 
in  this  area,  our  work  herein  exhibits  the  following  novelties:  1)  different  from 
[12,  13,  14],  which  consider  fixed  spreading  code  assignment,  we  introduce  one 
additional  degree  of  freedom,  i.e. ,  we  attempt  to  optimize  the  spreading  codes 
at  the  physical  layer  to  further  improve  the  system  performance.  We  develop 
a  novel  low-complexity  suboptimal  solution  to  an  otherwise  NP-hard  problem 
with  excellent  performance  in  practice;  2)  in  our  work,  interference  between  sec¬ 
ondary  links  is  taken  into  consideration,  and  the  SINR  requirement  necessary 
to  guarantee  a  certain  level  of  quality  of  service  for  secondary  users  is  also  guar¬ 
anteed.  This  is  different  from  [15]  and  [16],  where  only  the  interference  from 
secondary  to  primary  users  is  considered;  3)  direct  sequence  spread  spectrum 
technology  is  adopted  in  our  work,  while  [17]  and  [22]  consider  CDMA-based 
multi-carrier  modulation;  4)  no  previous  work  [12]  -  [22]  considers  interactions 
between  the  higher  layer  functionalities  of  the  networking  protocol  stack  (i.e., 
routing  and  scheduling)  and  physical  layer  resource  allocation  (i.e.,  power  and 
spreading  code  design)  in  a  multihop  cognitive  radio  network  scenario.  In  our 
work,  we  explore  a  new  framework  that  captures  the  interdependencies  between 
spread-spectrum  channelization  and  routing  in  cognitive  radio  ad  hoc  networks. 
To  the  best  of  our  knowledge,  our  work  is  the  first  to  consider  joint  code-channel 
optimization  and  routing  in  cognitive  ad  hoc  networks. 

C.  System  Model 

We  consider  a  primary  spread-spectrum  system  with  processing  gain  (code  se¬ 
quence  length)  G.  Denote  VIA  as  the  set  of  active  primary  communication  links. 
Let  (l,  k )  denote  the  link  with  transmitter  l  and  receiver  k.  Note  that  link  (l,  k) 
is  distinct  from  link  ( k,l ).  Each  primary  link  is  pre-assigned  with  an  unique 
code  sequence,  i.e.,  s for  link  (k,l).  We  let  A/"  =  {1,  -  ■  ■  ,N}  represent  a  fi¬ 
nite  set  of  secondary  users  (also  referred  to  as  nodes).  Secondary  users  do  not 
have  any  pre-assigned  code  sequence  and  opportunistically  send  their  data  by 
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optimizing  code  sequence  and  power. 

Traffic  flows  of  secondary  users  are,  in  general,  carried  over  multi-hop  routes. 
Let  the  traffic  demands  consist  of  a  set  T>  =  {1,  2,  •  •  •  ,  D}  of  unicast  sessions. 
Each  session  d  £  T>  is  characterized  by  the  destination  node  d,  d  £  Af  for  the 
traffic.  We  indicate  the  arrival  rate  of  session  d  at  node  i  as  nf(t)  at  time  t.  Each 
node  maintains  a  separate  queue  for  each  session  d  for  which  it  is  either  a  source 
or  an  intermediate  relay.  At  time  slot  t,  define  Qf(t)  as  the  number  of  queued 
packets  of  session  d  waiting  for  transmission  at  secondary  user  i.  Define  rf3  (t) 
(in  packets/s)  as  the  transmission  rate  on  fink  (i,j)  for  session  d  during  time 
slot  t,  and  R  as  the  vector  of  rates.  Note  that  nf(t)  represents  the  exogenous 
traffic  arrivals  at  node  i,  while  k=£i  rki(t )  represents  the  total  endogenous 

traffic  arrivals  at  node  i  resulting  from  routing  and  transmission  decisions  from 
other  nodes  k. 

For  V*  £  J\f ,  the  queue  is  updated  as  follows: 


Qf(t+  !) 


q?w  +  E  4w-E 


D.  Physical  Layer  Model 

We  recall  that  the  secondary  cognitive  ad  hoc  network  coexists  with  the  primary 
system  over  the  primary  licensed  band.  In  general,  the  transmitted  spread 
spectrum  signal  for  the  link  (i.  j)  is  denoted  by 

OO 

Uij(t)  =  E  \44dj  (t  -  mT )  ej(27r/°t+^)6i(m)  (1) 

771=1 

where  bi(m )  £  {— 1,+1}  is  the  mth  data  bit  (binary  phase-shift-keying  data 
modulation),  Ei  is  the  total  transmission  energy,  and  cj>i  is  the  carrier  phase 
with  carrier  frequency  fc;  s-ij {i)  is  the  normalized  unit-energy  user  signature 
waveform  with  duration  T  given  by 


G-i 

sij{t)  =  'y  ^  ~  3^c)  (2) 

3=0 

where  s ij(g),  g  =  0, 1, . . . ,  G  —  1,  is  the  value  of  the  gth  chip  of  the  spreading- 
code  vector  of  the  link  ( i,j ),  ip(t)  is  the  chip  waveform,  and  Tc  =  is  the  chip 
period.  The  combined  received  signal  waveform  due  to  all  link  transmissions 
over  flat  fading  channels  of  impulse  response  hij  (t) 

V(t)=  E  -  Tij)  +  n{t)  (3) 

( i,j)£VUuSU 

where  SU  is  the  set  of  active  secondary  finks,  £  [0,T)  is  the  relative  delay  of 
fink  (i,j)  with  respect  to  the  link  of  interest,  and  n(t)  is  a  white  Gaussian  noise 

8 


Approved  for  Public  Release;  Distribution  Unlimited 


process.  After  carrier  demodulation,  chip-matched  filtering  and  sampling  at  the 
chip  rate  over  the  duration  of  a  symbol  (bit)  period  of  G  chips,  the  received 
signal  at  primary  node  k  over  the  link  of  interest  (/,  k),  denoted  as  yjy ,  can  be 
represented  as 


y  ik 


Ik  — \/ ElLlkSlkbl  -f-  ^  (  y/ EmLmkSmnbn 


(m,i>)era 

+  ^  ^  y/ Eu  Luk  Cuv  bu  “h  life* 

(u,v)£SVi 


(4) 


We  note  that  the  superscript  “p”  in  (4)  denotes  primary  node.  Similarly,  let 
y\^)  denote  the  secondary  signal  received  at  node  j  over  the  link  of  interest  (i,  j) 


y\f  =y/EzLZJ  Cijbi  +  E  V  EmL  mj  S?i 

(m,n)^VU 

+  E  y/ Eu  Luj  C'u.'i:  bu  A  U  /  , 


(5) 


where  Ei  >  0,  bi  £  {±1},  and  s ^  £  RG,  ||sifc||  =  1,  denote  bit  energy,  informa¬ 
tion  bit,  and  normalized  signature  vector  of  primary  user  l  over  primary  link 
(l,k),  ( l,k )  £  VU ,  respectively;  Ei  >  0,  bi  €  {±1},  and  c ^  £  RG,  ||c,j||  =  1, 
denote  the  bit  energy,  information  bit,  and  normalized  signature  vector,  respec¬ 
tively,  of  secondary  user  i  over  secondary  link  ( i,j )  £  SU;  Lik ,  ( l ,  k)  £  VU 

and  Lij ,  (i,j)  £  SU  denote  the  path  coefficients  for  primary  link  ( l ,  k)  and  sec¬ 
ondary  link  (i,j),  respectively,  n^.  and  n;  represent  additive  white  Gaussian 
noise  (AWGN)  at  primary  node  k  and  secondary  node  j,  correspondingly,  inde¬ 
pendent  from  each  other  with  0  mean  and  autocovariance  matrix  a2 1.  We  note 
that  the  superscript  “s”  in  (5)  denotes  secondary  node. 

Rather  than  being  fixed,  the  interference  between  the  secondary  and  pri¬ 
mary  links  varies  with  the  transmit  power  and  spreading  code  allocation.  As 
modeled  in  (4)  and  (5),  interference  from  both  primary  links  and  secondary 
links  is  considered  in  our  design.  Specifically,  as  shown  in  (4),  which  is  the 
received  signal  at  the  receiver  node  k  of  primary  link  ( l,k ),  the  second  term 
(m,n)e™  \/ EmLmkSmnbm  represents  the  interference  from  other  primary 

links,  and  the  third  term  v\eSU  \J EuLukCuvbu  represents  the  interference 
from  secondary  links.  Similarly,  in  (5),  which  is  the  received  signal  at  receiver 
node  j  of  secondary  link  (i,  j),  the  second  term  J2(m,n)&vu  VErnLmjSmnbm  rep¬ 
resents  the  interference  from  primary  links,  and  the  third  term  (U,v)esu  y/EuLUj 

(“>  v)£{i,j) 

represents  the  interference  from  other  secondary  links. 

The  linear  filters  at  the  primary  and  secondary  receivers  that  exhibit  maxi- 
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mum  output  SINR  [24]  can  be  found  to  be 

wmaxsiNR,ifc  =  ciA^  sik,  (l,k)  G  VU , 

^maa;SINR,ij  —  Cy  .  ( / ,  j )  G  SU , 

where  Aifc  =  ElylfcME}’  Ay  =  E{y^s)yy)T},  ci  and  c2  are  arbitrary  positive 
constants,  i.e. ,  ci,C2  >  0.  Maximum  output  SINRs  at  the  receiver  end  for  link 
(l,k)  and  (i,j)  are,  respectively,  given  by 

SINRifc  =  L2lkEi sfkA^k)Blk,  V(Z,  k)  G  VU, 

SINR,;,  =  /.^/-',s/;Ax;;  y  s,r  V(f,  j)  G  5W, 

where  and  Ay^j)  are  the  disturbance  autocorrelation  matrices  at  the 

receiver  end  for  link  (7,  fc)  and  (i, ;/),  respectively,  defined  by 


A\0,fc) 


A 


E 


^-'mk^msmnsmn 


E 


LukE u 


( m,n)e'PU  (u,v)€lSIA 


A \(*,j) 


A 


^  '  LmjEmsmnSmn 


T  2  -rp  T  I  2t 

*-*11  j *-Ju'-’UV'-’uv  '  v 


{m,n)€.VU 


(«,«)esw 


C-  ■ 

In  our  cognitive  radio  setup,  the  normalized  channel  capacity  -^7  of  secondary 
link  (i,  j)  as  a  function  of  SINRy  is  given  by  [25] 

^  =  log2(l  +  ^5/AJRy)  (6) 

where  Cy  is  the  channel  capacity  and  kb  is  the  bandwidth  of  the  primary 
licensed  band.  Given  the  fixed  bandwidth  W,  the  channel  capacity  Gy  can  be 
evaluated  by  (6)  for  any  instantaneous  value  of  SINRy. 

III.  Methods,  Assumptions,  and  Procedures 

A.  Problem  Formulation 

Our  goal  is  to  design  a  distributed  cross-layer  control  algorithm  to  maximize 
the  secondary  network  throughput  by  jointly  and  dynamically  allocating  routes, 
code  sequence  and  transmit  power  for  each  secondary  link  along  the  path.  De¬ 
note  SU(t)  as  the  set  of  secondary  links  chosen  for  activation  during  time  slot  t, 
c(t)  =  {c ij(t)  :  ( i,j )  G  SU{t)}  and  E(f)  =  {E,(t)  :  (i,j)  G  SU{t)}  as  the  sets  of 
code  sequences  and  power  allocation  decisions  for  every  active  secondary  link. 
An  ideal  throughput-optimal  network  controller  should,  at  each  decision  period 
(i.e.  time  slot  t),  find  SU(t),  c(f),  and  E(Z)  to  maximize 

E  E  Cy(c(i),E(f))-AQy(f),  (7) 

i£SU  jtzSIA ,j^i 
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where  A Qij{t)  =  maxdex>  [Qf(t)  -  Qj(t)] 
given  by 

E(t)) 

0 

where  d^(t)  =  argmaxdeC  {Qf(t)  -  Qlj(t)}. 

The  objective  function  (7)  is  defined  based  on  the  principle  of  dynamic  back¬ 
pressure,  first  introduced  in  [26].  It  can  be  proven  [27]  that  a  control  strategy 
that  jointly  assigns  resources  at  the  physical/link  layers  and  routes  to  maximize 
the  weighted  sum  of  differential  backlogs  (with  weights  given  by  the  achievable 
data  rates  on  the  link)  achieves  throughput  optimality,  in  the  sense  that  it  is 
able  to  keep  all  network  queues  finite  for  any  level  of  offered  traffic  within  the 
network  capacity  region.  Moreover,  a  desirable  solution  should  enable  secondary 
users  to  dynamically  utilize  the  available  spectrum  resource  in  the  code  domain 
to  provide  SINR  guarantees  for  both  primary  and  secondary  users.  The  problem 
can  be  expressed  as 

PI  :  Find  :  SU{t),  E(t),  c(f) 

Maximize  :  E  E  Cij(c(t),E(t))  ■  A Qij(t) 

i£SU  jtzSU 

Subject  to  : 

SINRZfe  >  SINRpy,  V(Z,  k)  G  VU, 

SINR„  >  SINRcJp,  V(i,  j)  G  SU(t), 
ci:j(t)T Cij{t)  =  1,  Ei  <  Emax,  V(i,j)  G  SU{t). 

Therefore,  ideally,  a  throughput-optimal  policy  should  continuously  (i.e. ,  at 
each  time  slot)  assign  resources  on  each  network  link  by  solving  problem  PI  to 
optimality.  However,  exact  solution  of  PI  requires  global  knowledge  of  all  fea¬ 
sible  allocations  and  a  centralized  algorithm  to  solve  a  mixed  integer  non-linear 
problem  (NP-hard  in  general)  such  as  PI  on  a  time-slot  basis.  This  is  clearly  not 
practical  for  real-time  decision  making.  This  provides  the  rationale  for  our  dis¬ 
tributed  algorithm,  which  is  designed  to  provide  an  approximate  solution  to  PI 
based  on  real-time  distributed  decisions  driven  by  locally  collected  information. 
Note  that,  for  the  sake  of  simplicity,  we  will  drop  all  time  dependencies. 

In  the  following  sections,  we  first  discuss  the  decentralized  solution  for  joint 
routing  and  code-division  channelization  in  Section  V.  Specifically,  we  described 
in  detail  about  how  next  hops  are  selected  together  with  physical  layer  resource 
optimization  in  our  proposed  algorithm.  We  also  discussed  how  nodes  learn 
about  the  environment  to  make  distributed  decisions  on  routing  based  on  a  com¬ 
bination  of  physical  sensing  and  of  local  exchange  of  information  through  control 
packets  at  MAC  layer.  Then,  we  present  in  detail  our  power  and  spreading  code 
allocation  algorithm  to  maximize  the  SINR  for  secondary  links  in  Section  VI. 


.  The  transmission  rates  are  then 

if  d  =  d*j{t) 
otherwise 
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B.  ROCH:  Distributed  Joint  Routing  and  Code-Division 
Channelization 


We  now  present  the  decentralized  joint  routing  and  code-division  channeliza¬ 
tion  solution,  which  aims  at  maximizing  throughput  through  joint  opportunistic 
routing,  dynamic  code  sequence  optimization  and  transmit  power  control.  In 
the  proposed  solution,  backlogged  nodes  i  first  maximize  their  local  objective 
function  Ctj  ■  AQij  over  all  feasible  next  hops  j  by  optimizing  c ^  and  Ei  based 
on  locally  collected  spectrum  information  -  details  are  given  in  what  follows. 
Then,  in  case  of  channel  access  contention,  each  node  will  access  the  channel 
with  a  probability  that  is  a  monotonically  increasing  function  of  its  local  utility. 
We  now  describe  the  details  of  the  proposed  solution.  Every  backlogged  node  i 
performs  the  following  algorithm: 

1.  Find  the  set  of  feasible  next  hops  {nf ,  nf,  •  •  •  ,  nf}  for  the  backlogged  ses¬ 
sion  d,  which  are  neighbors  with  positive  advance  towards  the  destination 
of  d.  We  say  node  j  has  positive  advance  with  respect  to  i  iff  j  is  closer  to 
the  destination  d  than  i. 

2.  For  each  candidate  next  hop  j  £  {nf ,  nf,  ■  ■  ■  ,  nf},  maximize  link  capacity 
Cij  by  optimizing  Cij  and  Ei ,  using  the  algorithm  proposed  in  the  following 
Section  VI. 


3.  Schedule  the  session  d*  with  next  hop  j*  with  maximal  Cij  ■  AQij.  Hence, 
routing  is  performed  in  such  a  way  that  lightly  backlogged  nodes  with 
higher  link  capacity  receive  most  of  the  traffic. 


4.  Once  the  next  hop  and  corresponding  code  sequence  and  power  allocation 
have  been  determined,  the  probability  of  accessing  the  medium  is  calcu¬ 
lated  based  on  the  value  of  C',;?  ■  AQij.  Nodes  with  higher  Cjj  ■  AQ^  will 
have  a  higher  probability  of  accessing  the  medium  and  transmit.  Note  that 
links  with  higher  differential  backlog  may  have  higher  spectrum  utility,  and 
thus  have  higher  probability  of  being  scheduled  for  transmission.  This  is 
implemented  by  varying  the  size  of  the  contention  window  at  the  MAC 
layer.  The  transmitting  node  i  generates  a  backoff  counter  BCi  chosen 
randomly  (with  a  uniform  distribution)  within  the  interval  [1,  2CWi~1], 
where  CWi  is  the  contention  window  of  transmitter  i,  whose  value  is  a 
decreasing  function  4>()  of  Cij  ■  AQij  as  below 


CWt  =  -0i 


Cij  *  AQij 

£fc6V,lWev(C«  ‘  ^-Qkl) 


+  02, 


(9) 


where  0i  >  0,02  >  0  and  the  denominator  represents  the  objective  value 
of  neighboring  competing  nodes.  Note  that  sender  i  collects  the  utility 
values  of  its  neighbors  by  overhearing  control  packets  coded  by  a  common 
spreading  code. 


Note  that  only  the  second  step  involves  solving  a  real  optimization  (maxi¬ 
mization)  problem  (with  low  complexity).  Specifically,  in  step  2,  the  user  needs 
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Algorithm  1  ROCH  Algorithm 
1:  At  backlogged  node  i 
2:  for  each  backlogged  session  cl  do 
3:  for  j  G  {nf,  n%, ...,  n^}  do 

4:  Calculate  link  capacity  Cij  by  optimizing  c,y  and  E3  using  algorithm  2 

proposed  in  the  following  Section  VI. 

5:  end  for 

6:  end  for 

7:  Schedule  ( d*,j *)  =  argma x(C^  •  A Qtj) 

8:  Set  contention  window  CWj  =  —  ^ - ctJ-AQ,3  _p  g 

1  1  EfceA/yMevlCV'AQfci)  ^ 

9:  Return  [s*  ,j* ,CW%] 


to  solve  the  optimal  power  and  spreading  code  allocation  problem  by  performing 
the  algorithm  discussed  in  Section  VI  (which  has  polynomial  complexity).  In 
step  3,  the  user  simply  selects  the  next  hop  (and  the  associated  session)  with 
maximal  utility  by  comparing  utility  values  of  all  candidate  next  hops.  In  step 
4,  the  user  simply  calculates  the  contention  window  size  based  on  (9)  for  MAC 
layer. 

The  algorithm  calculates  the  next  hop  j  opportunistically  depending  on 
queueing  and  spectrum  dynamics,  according  to  the  objective  function  Cij  ■  AQ,j . 
The  combination  of  opportunistically  selected  next  hops  leads  to  a  multi-hop 
path.  The  multi-hop  path  discovery  terminates  when  the  destination  is  selected 
as  the  next  hop.  If  the  destination  is  in  the  transmission  range  of  the  transmitter 
(either  a  source  or  an  intermediate  hop  for  that  session) ,  the  differential  backlog 
between  the  transmitter  and  the  destination  is  no  less  than  the  differential  back¬ 
logs  between  the  transmitter  and  any  other  nodes,  because  the  queue  length  of 
the  destination  is  zero.  With  this  scheme,  lightly-congested  nodes  (as  indicated 
by  a  smaller  differential  backlog)  have  a  higher  probability  of  being  selected 
as  intermediate  relays.  Links  with  larger  differential  backlogs  result  in  smaller 
contention  window  size  at  the  medium  access  control,  and  therefore  have  higher 
probability  of  accessing  the  channel  to  reserve  resources.  Ultimately,  heavily 
backlogged  nodes  with  potential  high-data  rate  opportunities  to  transmit  have 
a  higher  probability  of  accessing  the  channel. 

The  algorithm  is  implemented  through  a  MAC  protocol  that  uses  a  three- 
way  handshake  between  source  and  destination.  The  three-way  handshake  is 
carried  out  via  exchange  of  Request-to-Send  (RTS),  Clear-to-Send  (CTS)  and 
Data  Transmission  reservation  (DTS)  packets.  Backlogged  nodes  contend  for 
spectrum  access  on  a  common  control  channel  that  operates  in  parallel  to  the 
data  channel  through  a  common  spreading  code.  Each  node  makes  adaptive 
decisions  based  on  local  information  collected  through  RTS/CTS/DTS  packets 
coded  with  a  common  spreading  code.  The  DTS  packet  is  used  to  announce 
the  information  on  spreading  code,  transmit  power  allocation  and  queueing  in¬ 
formation  to  neighboring  nodes.  We  include  queuing  information  in  the  control 
packets  to  allow  the  nodes  to  make  routing  decisions.  The  set  of  feasible  next 
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hops  can  be  obtained  by  requiring  only  the  neighborhood  geographical  informa¬ 
tion  (can  be  obtain  by  GPS  for  example)  and  an  estimate  of  the  destination’s 
position.  Geographical  information,  together  with  the  queue  information,  are 
encapsulated  in  the  control  packets  (i.e. ,  RTS/CTS/DTS)  to  allow  nodes  to 
exchange  their  information. 

C.  Cognitive  Secondary  Power  and  Spreading  Code  Allo¬ 
cation 

We  now  come  to  the  spread-spectrum  channelization  problem,  which  is  a  crucial 
component  of  our  spectrum  management  framework  and  is  executed  at  step  2 
of  the  algorithm  described  in  the  previous  section.  We  would  like  to  dynami¬ 
cally  and  adaptively  generate  digital  waveforms  that  span  the  whole  available 
spectrum  band  to  maximize  the  capacity  of  the  link,  and  that  at  the  same  time 
avoid  interfering  with  existing  users  (what  we  refer  to  as  spread-spectrum  chan¬ 
nelization).  Each  secondary  user  optimizes  the  waveform  to  be  used  based  on 
the  current  interference  condition  at  the  secondary  receiver  and  on  the  interfer¬ 
ence  receiver  profile  to  maximize  the  overall  network  throughput,  as  described 
in  Section  V.  We  present  the  secondary  code-channel  optimization  formulation 
in  Section  VI-A.  The  formulated  problem  is  NP-hard  with  complexity  exponen¬ 
tial  in  G.  In  Section  VI-B,  we  develop  a  realizable  suboptimal  solution  with 
polynomial  complexity. 

D.  Cognitive  Spread-spectrum  Channelization  Formula¬ 
tion 

In  order  for  a  cognitive  radio  network  to  efficiently  share  the  licensed  band  with 
the  primary  network,  the  secondary  transmitter  has  to  guarantee  the  SINR  QoS 
of  all  primary  users  and  ongoing  secondary  users  and,  simultaneously,  maximize 
the  secondary  SINR  of  link  ( i,j )  at  its  receiver  end.  We  first  assume  that 
the  secondary  link  (i,  j)  is  activated  with  the  transmission  bit  energy  Ei  and 
(normalized)  signature  vector  c,3 .  In  this  spirit,  our  objective  is  to  find  the  pair 
( Ei,Cij )  that  maximizes  SINR^  under  the  constraints  that  SINRjfc,  V(Z,fc)  £ 
VU,  and  SINRU„,  V(u,  v)  £  SU  —  {{i,j)},  are  all  above  the  prescribed  thresholds 
SINRpy  and  SINR.^p ,  respectively,  i.e. 

(Ei,Cij)opt  =  argmax  EicJjAV]  ...  dj  (10) 

Bi>0,ci36KG 

s.t.  EiL2lk sJkA^lk)  slk  >  SINR*pV  V(Z,  k)  £  VU, 

EuL2uv cLAC^)  cuv  >  SINRsp,  V(u,  v)  £  SU  -  {(i,j)}, 

Cij^ij  =  1,  Ei  ^  Emax  5 

where  Emax  denotes  the  maximum  allowable  bit  energy  for  the  secondary  user. 

The  optimization  task  of  maximizing  a  quadratic  objective  function  (A^  ^ 
is  positive  definite)  subject  to  the  constraints  in  (10)  is,  unfortunately,  a  non- 
convex  NP-hard  optimization  problem  [29].  In  the  following,  we  delve  into  the 
details  of  the  problem  and  develop  a  realizable  solution. 
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E.  Cognitive  Secondary  Channel  Design 

Using  the  matrix  inversion  lemma  on  A ^  k ^  and  Ad  %,  respectively,  we  can 
express  the  key  quadratic  constraint  expressions  in  (10),  respectively,  as 


sikA\(i,k)  Slk  ~ 

T  A  r1 


sTkAiksik 


-*■  ElL ‘iksTkAlk  Slk 


mT  A  -1 


-1_ 

Ik 

C  a 


,  V(Z,  k)  g  vu , 


/\  r  —  _ ~UV*‘-UV  '-UV 

^uvA\(u,v)  1  F  T2  CT  A-U 


-Ju-LJuvK'uvjrs-'uv  ^UV 


V(u,v)  g  SU-  {(*,  j')}, 


(11) 


(12) 


where  we  recall  that  Ajfc  =  E{y^y[^T}  and  Auv  =  E{yisJyisJT}-  Then,  the 
primary  SINR  constraints  in  (10)  become 


*TkAiksik  > 


SINR  % 


A 


EiL2lk 


SINR  %EiLl 


=  alk ,  V(l,k)eVU, 


Ik 


(13) 


while  the  secondary  SINR  constraints  are  given  by 


T  i  ^  SINR&,  a 

C“1'  UvCuV  ~  EuL*v  +  SINR %EuLlv  ~  luv ’ 

V(u,w)  g  5W-  {(*,  j)}- 


The  optimization  problem  (10)  can  be  rewritten  as 

(Ei,Cij)opt=  argmax  Eicf.Ar,] ...  dj 

Bi>o,cijeiG 

s.t.  sTkATksik^  aiki  v(/,/c)  g?w, 

cuvAuvcuv  ^  V(u,  u)  g  iSTf  ~  {(bj)}> 

cTc  =1  E  <  E 


(14) 


(15) 


Using  the  matrix  inversion  lemma  on  A^,  and  A“^  this  time,  we  can  ex¬ 
press  the  optimization  constraints  as  explicit  functions  of  signature  vector  of 
the  secondary  link  ctj ,  i.e. 


sTkAlkSl‘k  ^ 


lk\(i,j)CiiCijAlk\(i,j)Slk 


1  +  EiL?kcf;jAlk^ijjCij 


OLlk, 


V(Z,  fc)  g 


(16) 


cT  A_1c  > 


-1 


1  +  EiL? ,ci,  A 


-l 


ikK'ijJr^u.v\(i,j)Cij 


V(u,u)  g  SW-  {(*,  j)}- 


(17) 
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For  notational  simplicity,  define  the  G  x  G  matrices 


where 


B,,  =  L*A-{ 


B  =  L2  A”1  c  cT  A-1  —  3  T?  A-1 


filk  —  ^lkA-lk\(i,j)Slk  ^Iki  V(z,  k)  G  PW, 


(18) 

(19) 

(20) 


/3UV  =  CuvA-uv\(itj)cuv  ~  luvi  y(u,v)  £  SU  —  {(*,j)}-  (21) 

Then,  the  optimization  problem  in  (15)  can  be  rewritten  -for  one  more  time-  as 

xopt  =  arg  max  x1  Ar,1  x  (22) 

X£RG  V  J 

s.t.  xTB;fcx  -  Pik  <  0,  V(Z,  k)  £  VU, 

xTBuvx  -  puv  <  0,  V(u,  v)  £  SU  -  {(i,j)}, 

X  X  Ernax 

where  x  is  the  amplitude-including  signature  vector  of  secondary  link  (i.  j), 

x  =  \J EjCjj .  From  the  perspective  of  computational  effort,  however,  (i)  B ik 
and  Buv  are  not  necessarily  positive  semidehnite,  hence  the  problem  in  (22) 
is  in  general  a  non-convex  quadratically  constrained  quadratic  program  (non- 
convex  QCQP),  and  (ii)  the  complexity  of  a  solver  of  (22)  is  exponential  in  the 
dimension  G  (NP-hard  problem). 

Using  the  commutative  property  of  trace,  we  are  able  to  represent  the  opti¬ 
mization  problem  in  (22)  as 

Xopt  =  arg  max  lf{A:‘  X}  (23) 

xeRGxc  u 

s.t.  Tr{BlkX}  <  fak,  V(Z,  k)  £  VU, 

7V{BU„X}  <  0UV,  \/(u,v)  GSU  -  {(i,j)}, 

Tr{X}  <  Emax,  X  ^  0,  ranfc(X)  =  1. 

where  X  =  xxT  and  X  ^  0  denotes  that  the  matrix  X  is  positive  semidefinite. 

To  effectively  attack  the  problem  anyway,  the  problem  in  (23)  is  relaxed  to 
a  semidefinite  program  by  dropping  the  rank  constraint,  i.e. 

X' =  arg  max  Tr{  A~\  ,.X}  (24) 

xsrgxg  u  ,JI 

s.t.  Tr{BifeX}  <  plkt  V(Z,  k )  £  VU, 

Tr{BuvX}  <  fiuv,  W(u,v)  £  SU  -  {( i,j )}, 

Tr{X}  <  Emax,  X  >z  0. 
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Then,  (24)  is  a  convex  problem  that  can  be  solved  using  semidefinite  program¬ 
ming  within  the  error  e  in  0(G4  log  1/e)  time.  The  solution  X"  returned  by 
semidefinite  programming  makes  the  objective  function  Tr{A/)4^X}  attain 
a  value  within  (TrjAw^X'}  —  e,  Tr{ Ab  .jX'}).  Of  course,  because  of  the 
constraint  relaxation  itself  the  objective  function  evaluated  at  the  optimum 
point  X'  in  (24)  is  just  an  upper  bound  on  the  value  of  the  objective  function 
evaluated  at  the  optimum  point  of  interest  Xopt  of  (23),  Tr{A/d  ^Xopt}  < 
Tr{A/(4  j)X'}.  Since  X'  is  not  available,  we  have  X"  with  Tr{A//  yX"}  £ 
(Tr{Ar/^ :)X'}  —  e,  Tr{A/)4  ^X'}).  So  far,  for  the  cognitive  design  of  a  code¬ 
division  secondary  link,  first,  for  the  given  primary  and  secondary  SINR-QoS 
thresholds  SINRpp  >  0  and  SINR^p  >  0,  we  test  whether  fak,  V(Z,fc)  £  VU 
and  Puv,V(u,v)  £  SU  —  {(*,  j)}  are  all  greater  than  zero.  If  this  is  not  true, 
then  the  SINR-QoS  constraints  cannot  be  satisfied  and  outright  no  secondary 
transmission  is  allowed.  Otherwise,  we  solve  the  problem  (24)  by  semidefi¬ 
nite  programming  with  the  return  X".  If  the  rank  of  X"  is  1  with  eigenvalue, 
eigenvector  pair  (Ai,vi),  then  we  already  have  our  secondary  link  design  with 
signature  cl:j  =  V!  and  transmission  bit  energy  £/  =  Ai .  If  the  rank  of  X"  is 
not  1,  further  work  is  needed  as  described  below. 

When  X'  of  (24)  (or  in  practice  X")  is  of  rank  2  or  more,  we  can  not  find  the 
direct  mapping  from  X'  of  (24)  to  xopt  in  (22).  Under  this  case,  we  may  switch 
the  search  for  an  optimal  vector  in  (22)  to  a  search  for  an  optimal  probability 
density  function  (pdf)  of  vectors  that  maximizes  the  average  objective  function 
subject  to  average  constraints,  i.e. 

/op4(x)  =  max  E{xTAr(1  .}x}  (25) 

/(x)  U 

s.t.  E{xTB;fcx}  <  /3ik,V(l,  k)  £  VU, 

E{xtBm„x}  <  puv,V(u,v)  £SU  -  {( i,j )}, 

E{xtx}  <  Emax 

where  /(x)  denotes  the  probability  density  function  of  x  and  E{-}  denotes  sta¬ 
tistical  expectation.  Using  the  commutative  property  between  trace  and  ex¬ 
pectation  operators,  the  pdf  optimization  problem  in  (25)  takes  the  equivalent 
form 


/^(x)  =  argmax  Tr{ AT*  E{xxT}}  (26) 

/(X)  u  ' 

s.t.  B;fcE{xxT}  <  /3ik,V{l,  k)  £  VU, 

B„„E{xxr}  <  puv,V(u,v)  eSU  -  {( i,j )}, 

E{xxt}  <  Emax. 

We  can  show  that  /opt(x)  is  in  fact  Gaussian  with  0  mean  and  covariance 
matrix  X',  fopt(x)  =  A/’(0,X').  With  X"  as  a  close  approximation  of  X',  we 
can  draw  now  a  sequence  of  samples  x1,x2,...,xp  from  A/”(0,X").  We  test 
all  of  them  for  “feasibility”  on  whether  all  constraints  of  (22)  are  satisfied  and 

17 


Approved  for  Public  Release;  Distribution  Unlimited 


among  the  feasible  vectors  (if  any)  we  choose  the  one,  say  x*'0),  with  maximum 
xrAr,1j  ?,x  objective  function  value.  We  could  have  suggested  at  this  time  a 

cognitive  secondary  link  design  with  y/EiCij  =  x^0).  Instead,  we  will  use  x^0) 
as  an  initialization  point  to  an  iterative  procedure  below  that  will  lead  to  a 
improved  link  design  vector.  First  we  express  as 

=  SE  ST  +  ct2I  (27) 


where  S  =  [{sik,  V(Z,  k)  £  VU}:  {cuv,  V(u,  v)  £  SU  —  denotes  the  matrix 

with  columns  the  signatures  of  the  primary  users  and  ongoing  secondary  users; 

E  =  diag ({£i£?j,V(Z,fc)  £  VU},{EuL2uj,M(u,v)  £  SU  -  {(i,j)}).  Using  the 
matrix  inversion  lemma, 


a -1 


=  iI-4s(u-1 


1 


STS)_1ST. 


Substitution  of  (28)  in  the  objective  function  of  (22)  leads  to 

S(bi) 


xtA,  }.  ,x  =  ^xTx - ^j-xTQx 


(28) 


(29) 


where  Q  =  S(E_1  +  ^STS)_1ST.  In  (29),  the  first  term  ^xTx  is  a  convex 
function  while  the  second  term  —  yrX2  Qx  is  a  concave  function  (which  implies 
that  yrX2  Q  x  is  convex).  Based  on  the  first-order  condition  of  convex  functions, 
we  have 

xTx  >  2x(°)tx  -  x(°>Tx(°>  (30) 

where  x^0)  denotes  an  initial  feasible  vector.  Then,  we  combine  (29)  and  (30)  and 
form  an  optimization  problem  that  maximizes  the  following  concave  function 

2  X(0)TX  _  1  XTQ  x  _  J_x(0)Tx(0)  (31) 

oA  er4  az 

that  leads  to  a  suboptimum  solution  for  our  original  problem  in  (22).  To  maxi¬ 
mize  (31)  in  view  of  our  constraints  in  (22),  we  restrict  all  non-convex  constraints 
into  convex  sets  (linearization).  In  particular,  we  consider  the  non-convex  con¬ 
straints 

xtB ifc  x  <  /3|fc,  V(Z,fc)€Z$,  (32) 

xtB„„  x  <  V(u,  v)  £  XW,  (33) 

where  Inc  and  Inc  denote  sets  of  all  pairs  (l,k)  £  VU  and  (u,v)  £  SU  — 
{ (z,  j ) }  for  which  the  quadratic  constraint  is  non-convex  function,  respectively. 
Then,  we  decompose  the  matrices  B and  B.U1;  into  its  positive  and  nega¬ 
tive  parts  that  are  positive  semidefinite,  i.e.  B n-  =  B^.  —  Bjj.,  where  B^.  = 
B ik  =  B uv  =  ~  B uvt  where  B+,  = 

LljA\(ij)CuvCuvA\(i,j)>  B uv  =  PuvLl jA^-).  Therefore,  the  original  non- 
convex  constraints  can  be  written  as 


xtB+  x-/3ik<  xtB-  x,  V(Z,  k)  £  XW,  (34) 
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xTB+  x  -  /3UV  <  xtBuv  x,  V(it,  v)  G  2$ ,  (35) 

where  both  sides  of  the  inequalities  are  convex  quadratic  functions.  Lineariza¬ 
tion  of  the  right-hand  side  of  (34)  around  the  vector  x^0)  leads  to 

xtB+  x  -  Afc  <  x(°)tB-  x<°)  +  2x(°)TB-fc  (x  -  xM), 

(36) 


X  -  , 


<  x(0)tB,T„x(0)+2x(0)tB: 


(x-x(0)), 


V(«,«)et> 


(37) 


In  (36)  and  (37),  the  right-hand  side  is  an  affine  lower  bound  on  the  origi¬ 
nal  function.  It  is  thus  implied  that  the  resulting  constraints  are  convex  and 
more  conservative  than  the  original  ones,  hence  the  feasible  set  of  the  linearized 
problem  is  a  convex  subset  of  the  original  feasible  set.  Thus,  by  linearizing 
the  concave  parts  of  all  constraints,  we  obtain  a  set  of  convex  constraints  that 
are  tighter  than  the  original  non-convex  ones.  Now,  the  original  optimization 
problem  takes  the  form 


x« 

s.t. 


:  arg  max 


r(o  )T 


x - jXTQx 


R0)Ty(0) 


kj  vj 

tB+  X  -  x^Bj,  (2x  -  x<°>)  <  Afc,  V(i,  k)  G  2$, 
TBlkx<plk, 

tB+  x  -  x(°)tB-  (2x  -  xf°>)  <  /3„„,  V(«,  v)  €  2« 
1  B„„x  <  fiuv,  V(u,v)  G  SU  -1$  -  {(*,  j)}, 

TX  <  Er, 


(38) 


The  problem  in  (38)  is  a  convex  QCQP  problem  and  can  be  solved  efficiently 
by  standard  convex  system  solvers  [30]  to  produce  a  new  feasible  vector  x^1^. 
The  objective  function  x1  ArvJ.  ^x  in  (22)  evaluated  at  x^1)  takes  a  value  that  is 

larger  than  or  equal  to  its  value  at  x^0-* .  Repeating  iteratively  the  linearization 
procedure,  we  can  obtain  a  sequence  of  feasible  vectors  x(°\  x^1),  x^2), . . . ,  x^ 
with  non-decreasing  values  of  the  objective  function  in  (22).  As  demonstrated  by 
experimental  results  in  [31],  it  is  observed  that  eight  or  nine  iterations  are  enough 
for  effective  convergence.  After  numerical  convergence,  the  secondary  link  (i,j) 
is  suggested  with  signature  c^-  =  x^/Hx^H  and  bit  energy  Ei  =  ||x^Ar)||2. 
If  the  secondary  transmission  bit  energy  and  signature  vector  returned  by  our 
algorithm  satisfy  the  constraint  E{  L'f^  cf.-  A^l  ^ ,  c,;;-  >  SINR^j,  the  secondary 
link  (i.  j)  with  Ei  and  c ij  is  considered  as  candidate  for  transmission.  Otherwise, 
transmission  over  the  path  (i.  j)  is  not  allowed. 

Our  proposed  scheme  is  summarized  in  Algorithm  2. 

The  complexity  of  the  optimization  in  the  physical  layer  (i.e.  optimizing 
the  secondary  transmission  bit  energy  and  signature  vector)  is  dominated  by 
the  complexity  of  solving  the  semidefinite  program  in  (24).  The  computation 
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Algorithm  2  Cognitive  Code-division  Channelization. 

if  (3ik  >  0 y(l,k)  £  VU  and  >  0,V(u,v)  £  SU  —  {(i,j)}  then 
Run  SDP  optimizer 
if  Rank(X")=l  then 

Obtain  (Ai,  vi)  from  eigen  decomposition  of  X" 

Assign  Ei  ■£-  Ai  and  c ij  <—  vi 

else 

Draw  xi,  X2, . .  .  ,  xp  from  Af(0,  X") 

if  Any  feasible  sample  xp,  p  £  {1,...,P}  that  satisfies  xjB;fexp  V(l,k)  £ 
VU  &XpB u„xp  V  (u,v)  £  SU  —  { (i,  j)}  then 

x(°)  feasible  xp  with  maximum  xjAr,^xp  value 

Run  iteratively  linearized  optimizer  with  convergence  point  x(iV) 

Assign  Ei  ■£-  ||xWf  and  c ij  ■£-  x(Jvi/||x(Jvi || 
else 

No  feasible  solution  by  assigning  Ei  ■£-  0  and  dj  <—  0 

Return 
end  if 
end  if 

if  EiLljcJjA^  ^dj  >  SINRfu  then 
Output  the  solution  Ei  and  c ij 
else 

No  feasible  solution  by  assigning  Ei  ■£-  0  and  c ij  <—  0 

end  if 
else 

No  feasible  solution  by  assigning  Ei  4—  0  and  cy  <—  0 

end  if 

Return  the  link  design  (Ei,  c ij)  as  candidate  for  transmission. 
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Fig.  2:  Secondary  receiver  SINR  as  a  function  of  the  iteration  step  of  the  lin¬ 
earized  optimizer  initialized  at  the  best  feasible  sample  out  of  P  =  20  drawings 
from  the  Af(0,  X")  pdf. 


complexity  of  the  algorithm  for  the  physical  layer  is  0(G4  log  1/e).  [31]  gives 
more  detailed  complexity  analysis  on  the  optimization  in  the  physical  layer.  For 
any  of  the  K  feasible  next  hops,  each  node  executes  ROCH,  which  clearly  has 
polynomial  complexity  in  the  number  of  0(G4  log  1/e).  Conversely,  the  central¬ 
ized  algorithm  of  this  family  of  schedule  and  routing  has  worst-case  exponential 
complexity  [28]. 

Here,  we  demonstrate  the  effective  convergence  of  proposed  algorithm  through 
simulation  studies.  In  the  simulation,  the  transmission  SNRs  of  primary  links 
are  all  set  equal  to  15  dB.  The  number  of  active  primary  users  is  set  to  16 
(fully  loaded).  All  the  signatures  for  primary  links  are  generated  from  a  mini¬ 
mum  total-squared  correlation  optimal  binary  signature  set  which  achieves  the 
Karystinos-Pados  (KP)  bound  [32]-[34] .  The  SINR  threshold  is  set  to  3  dB  for 
both  primary  and  secondary  users.  In  Fig.  2,  we  plot  the  secondary  receiver 
average  SINR  for  the  experimental  instants  of  ranfc(X")  >  1  as  a  function  of 
the  iteration  of  the  optimizer  initialized  at  the  point/design  x^0)  that  is  the 
best  out  of  P  =  20  samples  drawn  from  the  A/”(0,X")  pdf.  It  is  pleasing  to 
observe  that  eight  or  nine  iterations  are  enough  for  effective  convergence  and 
the  improvement  is  more  than  5  dB. 
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IV.  Results  and  Discussion 


In  this  section,  we  evaluate  the  performance  of  ROCH  in  a  multi-hop  cognitive 
radio  network.  We  have  developed  an  object-oriented  packet-level  discrete-event 
simulator,  which  models  in  detail  all  layers  of  the  communication  protocol  stack, 
including  routing,  medium  access  control  and  spread-spectrum  channelization 
as  described  in  this  project. 

A.  Performance  of  Secondary  Dynamic  Spectrum  Access 

We  consider  K  active  links  (including  primary  and  secondary  links)  with  signa¬ 
ture  length  (system  processing  gain)  G  =  16.  We  are  interested  in  establishing 
a  secondary  code-division  transmitter /receiver  pair  when  K  varies  from  14  to 
20.  All  link  signatures  are  generated  from  a  minimum  total-squared-correlation 
optimal  binary  signature  set  which  achieves  the  Karystinos-Pados  (KP)  bound 
for  each  ( K,G )  pair  of  values3.  The  transmission  SNRs  of  K  active  links  are 
all  set  equal  to  ^  =  15  dB ,  i  =  1,  2, ...  K\  the  maximum  allowable  transmis¬ 
sion  SNR  for  the  secondary  link  is  set  to  Emxx  =15  dB.  The  channels  are 
modeled  as  Rayleigh  fading.  The  channel  coefficients  are  taken  to  be  the  mag¬ 
nitude  of  independent  complex  Gaussian  random  variables  with  mean  0  and 
variance  1.  The  receiver  SINR  thresholds  for  primary  and  secondary  links  are 
set  to  SINRpy  =  3  dB  and  SINR^  =  3  dB,  respectively.  When  random  vector 
drawing  is  necessary,  P  =  20  test  vector  points  are  generated. 

In  Fig.  3(a),  we  plot,  as  a  function  of  the  number  of  active  links  K,  the 
percentage  of  times  that  one  more  secondary  link  is  enabled  directly  under  the 
proposed  algorithm  in  Section  V  as  well  as  the  random  code  assignment  (RCA) 
scheme  where  the  candidate  signature  is  randomly  generated  with  unit  norm 
and  transmission  bit  energy  is  set  to  be  the  maximum  allowable  value  Emax. 
We  observe  that  our  proposed  algorithm  in  Section  VI  offers  more  opportunities 
for  cognitive  secondary  transmission  than  RCA. 

In  Fig.  3(b),  to  gain  visual  insight  into  operation  of  the  network  we  plot 
the  instantaneous  receiver  SINR  of  a  primary  active  link  and  the  candidate 
secondary  link  for  the  case  K  =  17  over  an  experimental  data  record  sequence 
of  1000  Rayleigh  fading  channel  realizations.  Missing  secondary  signal  SINR 
values  indicate  instances  when  no  feasible  solution  was  returned.  The  proposed 
scheme  highly  increases  the  probability  of  secondary  transmission  compared  to 
RCA.  When  secondary  transmissions  do  occur  for  both  schemes,  joint  power 
and  sequence  optimization  executed  by  the  proposed  scheme  results  in  superior 
SINR  performance  for  the  secondary  receiver  over  RCA. 

B.  Network  Performance  Evaluation 

A  grid  topology  of  49  secondary  nodes  and  14  active  primary  links  is  deployed  in 
a  5000  m  x  5000  m  area.  The  spreading  code  length  is  set  to  16  for  both  primary 

3 For  G  =  16,  when  K  <  G  the  KP-optimal  sequences  coincide  with  the  familiar  Walsh- 
Hadamard  signature  codes. 
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and  secondary  users.  Active  primary  links  use  pre-assigned  code  sequences  as 
described  in  the  previous  section.  We  initiate  CBR  traffic  sessions  between 
randomly  selected  but  disjoint  source-destination  pairs  among  the  49  nodes. 
Parameters  Q\  and  62  in  (9)  are  set  to  5  and  10,  respectively.  Rayleigh  fading 
channel  is  used  and  the  path  loss  exponent  is  set  to  four.  We  compare  the 
performance  of  ROCH  with  three  alternative  schemes.  All  alternative  schemes 
rely  on  the  same  knowledge  of  the  environment  as  ROCH.  In  particular,  we 
consider  the  solution  SP-SIG  where  routing  is  based  on  the  shortest  path  with 
dynamic  signature  and  power  allocation  (as  proposed  in  Section  VI) .  The  other 
two  solutions  use  RCA  with  fixed  power  allocation.  We  consider  BP-RCA  as 
the  solution  where  routing  is  based  on  the  same  utility  function  as  ROCH  but 
with  random  code  assignment,  and  SP-RCA  as  the  solution  where  routing  is 
based  on  the  shortest  path  with  random  code  assignment. 

We  compare  the  four  solutions  by  varying  the  number  of  sessions  injected 
into  the  network  and  plot  the  network  throughput  (sum  of  individual  session 
throughputs).  Figures  3(c)  and  4(a)  show  the  impact  of  the  number  of  sessions 
injected  into  the  network  on  the  network  throughput.  The  traffic  load  per  ses¬ 
sion  is  set  to  2  Mbit/s  and  4  Mbit/s.  As  shown  in  both  figures,  ROCH  achieves 
the  highest  throughput.  The  improvement  obtained  by  ROCH  is  more  visible 
when  the  number  of  active  sessions  increases.  Solutions  with  adaptive  signature 
design,  i.e.,  the  proposed  solution  ROCH  and  SP-SIG,  outperform  RCA-based 
solutions,  i.e.,  BP-RCA  and  SP-RCA.  With  the  same  signature  optimization  al¬ 
gorithm,  ROCH  outperforms  SP-SIG  since  SP-SIG  restricts  packets  forwarding 
to  the  receiver  that  is  the  closest  to  the  destination,  even  if  the  link  capacity  is 
very  low  or  the  receiver  is  heavily  congested. 

Fig.  4(b)  shows  the  impact  of  source  data  rate  per  session  on  the  performance 
of  throughput.  We  evaluate  the  throughput  performance  as  the  traffic  load  per 
session  increases  from  1  Mbit/s  to  4  Mbit/s.  As  shown  in  Fig.  4(b),  ROCH 
obtains  a  significant  throughput  gain. 

Fig.  4(c)  shows  the  delay  performance  for  the  four  solutions  with  traffic 
load  2  Mbit/s  per  session.  In  general,  solutions  with  adaptive  signature  de¬ 
sign  (ROCH  and  SP-SIG)  outperforms  RCA-based  solutions  (BP-RCA  and  SP- 
RCA)  in  terms  of  delay,  and  the  delay  performance  gap  between  the  two  grows 
as  the  number  of  sessions  increases. 
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Number  of  Primary  Users  K 


(a) 


Instantaneous  Receiver  SINR 


(b) 


Throughput  vs  #  of  Active  Sessions 
Source  Rate  per  Session  =  2  Mbit/s 


(c) 

Fig.  3:  (a)  Secondary  transmission  percentage  as  a  function  of  the  number  of 
active  links  under  Cases  rank(X.")  =  1  and  >  1  (the  study  includes  also  the 
random  code  assignment  scheme) ;  (b)  Instantaneous  output  SINR  of  a  primary 
signal  against  primary  SINR-QoS  threshold  SINRpp  (thick  line)  and  instanta¬ 
neous  output  SINR  of  secondary  signal;  (c)  Throughput  vs.  Number  Active 
Sessions. 
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Throughput  vs  #  of  Active  Sessions 


#  of  Active  Sessions 


(a) 


Throughput  vs  Traffic  Load 
8  Active  Sessions 


(b) 


Delay  vs  #  of  Active  Sessions, 


(c) 

Fig.  4:  (a)  Throughput  vs.  Number  Active  Sessions;  (b)  Throughput  vs.  Traffic 
Load;  (c)  Delay  vs.  Number  Active  Sessions. 
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C.  RcUBe:  Real-Time  Reconfigurable  Radio  Framework 
with  Self-Optimization  Capabilities 

As  of  today,  existing  wireless  networking  standards  are  inherently  hardware- 
based.  The  majority  of  deployed  networks  are  neither  capable  of  keeping  up  with 
the  quickly  evolving  pace  of  wireless  standards,  nor  of  adopting  application- 
specific  protocol  designs,  mainly  due  to  their  limited  programmability  and  re¬ 
configurability  [35,  36].  As  a  consequence,  practical  deployment  of  new  protocol 
amendments  or  of  new  standards  and  their  application  is  prohibitive  in  terms 
of  both  cost  and  time.  There  is  a  profound  and  urgent  need  for  reconfigurable 
wireless  architectures  to  overcome  the  aforementioned  constraints  and  offer  the 
required  flexibility. 

Moreover,  while  in  recent  years  there  have  been  significant  algorithmic  de¬ 
velopments  in  cross-layer  network  adaptation  and  resource  allocation  [37,  38,  39, 
40],  in  most  cases  existing  architectures  are  not  flexible  enough  to  incorporate 
these  advancements  other  than  through  ad  hoc  tweaks.  As  a  consequence,  algo¬ 
rithmic  developments  and  architectural  and  protocol  innovation  follow  parallel 
paths,  with  little  or  no  cross-fertilization  between  the  two  communities. 

For  example,  among  the  emerging  wireless  technologies,  cognitive  radio  (CR) 
has  been  envisioned  as  a  technology  with  a  strong  potential  to  increase  the 
spectrum  utilization  efficiency  [41].  Cognitive  radio  technology  relies  on  the 
ability  of  the  radio  terminal  to  reconfigure  crucial  communication  parame¬ 
ters  (e.g.,  modulation,  carrier  frequency,  transmission  power,  among  others) 
to  achieve  dynamic  spectrum  access.  In  recent  work  [42,  43,  44,  45,  46,  47, 
48,  49],  the  importance  and  need  for  reconfigurability  between  different  layers 
of  the  protocol  stack  (cross-layer  approach)  has  been  highlighted  and  demon¬ 
strated.  However,  it  is  unclear  how  these  developments  can  be  incorporated 
in  existing  and  future  network  architectures.  Moreover,  as  mentioned  in  [50], 
application-specific  (e.g.,  QoS-featured  applications  and  white-space  network¬ 
ing)  and  context-specific  (e.g.,  video  communications)  protocol  designs  require 
the  establishment  of  more  flexible,  reconfigurable  wireless  infrastructures.  Ev¬ 
idently,  reconfigurability  is  becoming  a  crucial  feature  for  wireless  networking 
technologies. 

Based  on  these  premises,  in  this  project  we  introduce  Real-time  Reconfigurable 
Radio  (RcUBe),  a  novel  architectural  radio  framework  that  offers  self-optimization 
and  real-time  reconfiguration  and  adaptation  capabilities  at  the  PHY,  MAC, 
and  network  layers  of  the  networking  protocol  stack.  RcUBe  introduces  an  ex¬ 
tensive  and  evolving  abstraction  for  implementing  a  wide  range  of  protocols. 
For  this  purpose,  RcUBe  adopts  the  idea  of  decomposing  protocols  into  prim¬ 
itive  building  blocks  [35,  36,  51,  50].  The  novelty  of  RcUBe  is  based  on  its 
extent  and  capabilities.  First,  unlike  [35,  36,  51,  50]  RcUBe  offers  real-time 
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reconfigurability  at  three  layers  (PHY,  MAC,  routing).  Second,  it  provides  a 
structured  methodology  to  implement  cross-layer  algorithms  and  optimizations 
spanning  multiple  layers  of  the  protocol  stack.  Finally,  it  provides  support  for 
the  implementation  of  real-time  self-optimization  capabilities  that  are  the  basis 
of  state-of-the-art  CR  algorithms. 

The  design  of  RcUBe  is  modular,  in  the  sense  that  it  preserves  a  layered  pro¬ 
tocol  stack;  and  is  compatible  with  TCP/IP  (although  integration  with  TCP/IP 
is  out  of  the  scope  of  this  project).  However,  the  architecture  of  the  network 
node  is  divided  into  four  distinct,  separate,  but  interacting  planes,  each  logi¬ 
cally  in  charge  of  a  different  group  of  functionalities.  These  planes  are  referred 
to  as  decision ,  control ,  data ,  and  register  planes.  The  decision  plane  consists 
of  user-defined  decision  algorithms  that  provide  the  control  logic  for  real-time 
protocol  optimization,  adaptation,  and  reconfiguration.  In  the  control  plane, 
RcUBe  stores  the  required  logic  for  routing  and  MAC  protocol  execution  as 
well  as  data  plane  management.  Data  processing  takes  place  in  the  data  plane, 
while  the  register  plane  stores  and  manages  access  to  system  parameters  and 
environmental  information.  This  division  has  several  advantages.  First ,  this 
structure  separates  the  decision  plane,  which  enables  real-time  optimization, 
from  execution  of  the  protocol  stack.  Thus,  it  becomes  possible  to  define  and 
reconfigure  the  decision  logic  on-the-fly  without  influencing  the  on-going  pro¬ 
tocol  execution  logic.  The  decision  plane  consists  of  a  decision  engine  where 
user-defined  algorithms  are  executed  to  enable  adaptation  and  optimization  of 
the  protocol  stack.  Second ,  we  decouple  data-processing  functionalities  from 
their  control  logic.  In  other  words,  we  decompose  the  protocol  stack  into  data 
and  control  planes.  In  this  way,  we  obtain  the  essential  modularity  and  flex¬ 
ibility  for  reconfiguring  the  protocol  stack  in  real-time.  Finally ,  we  define  a 
register  plane,  which  represents  a  shared,  common  storage  space  of  the  state 
variables  of  protocols  handling  functionalities  at  different  layers  of  the  protocol 
stack.  The  register  plane  provides  a  liaison  between  different  planes  and  also 
stores  system  parameters  and  general  purpose  variables  that  can  be  accessed  by 
multiple  protocols. 

Our  design  is  implemented  and  tested  on  the  Universal  Software  Radio  Pe¬ 
ripheral  (URSP)  platform  [52],  a  commercial,  low-cost  RF  front-end.  This  is 
in  contrast  to  [51]  that  uses  comparatively  expensive  WARP  [53]  boards.  C+- 1- 
and  Python  languages,  supported  by  the  GNU  Radio  framework,  are  used  for 
PHY  layer  implementations,  while  the  overall  architecture  is  implemented  in 
Python.  Through  these  design  choices,  RcUBe  offers  easier  programmability 
compared  to  [35,  51,  36].  The  proposed  framework  is  natively  decoupled  from 
hardware.  In  its  current  instantiation,  it  mainly  runs  on  a  host  PC,  but  it 
can  be  adapted  to  run  on  any  multi-core  platform  along  with  any  available  RF 
front-end;  as  well  as  on  soft-core  processors  implemented  on  FPGAs.  RcUBe  is 
based  on  a  software-defined  physical  layer,  so  it  provides  broad  reconfigurabil¬ 
ity  at  the  PHY  layer,  as  it  is  not  strictly  limited  to  the  reconfiguration  of  the 
capabilities  of  selected  commodity  WLAN  cards.  Instead,  it  exploits  the  PHY 
layer  capabilities  of  the  GNU  Radio  framework  and  leverages  its  modular  and 
flexible  structure. 


27 


Approved  for  Public  Release;  Distribution  Unlimited 


To  sum  up,  the  major  contributions  of  this  project  are: 

•  the  design  of  a  real-time  reconfigurable  radio  framework  with  support  for 
control  logic  enabling  adaptation  and  self-optimization; 

•  real-time  reconfigurability  for  PHY,  MAC,  and 
routing  protocols; 

•  real-time  cross-layer  optimization; 

•  an  extensive  abstraction  for  implementing  a  wide  range  of  protocols; 

•  illustration  of  the  unique  features  and  the  flexible  design  of  RcUBe  by 
implementing  a  cross-layer  architecture  that  performs  on-the-fly  recon¬ 
figurations  based  upon  decisions  taken  in  real  time  by  a  decentralized 
optimization  algorithm. 

D.  Related  Work 

Wireless  protocol  stacks  are  traditionally  implemented  on  hardware  with  a  static 
architecture,  which  offers  little  or  no  room  for  reconfigurability.  To  overcome 
this  limitation  and  enable  flexibility  on  traditional  wireless  infrastructure,  early 
work  adopted  the  idea  of  using  multiple  existing  modules  on  commercial  802.11 
cards.  For  example,  MultiMAC  [54]  is  based  on  the  idea  of  switching  between 
pre-defined  MAC  protocols.  However,  this  offers  only  limited  flexibility,  since 
only  a  limited  number  of  MAC  protocols  can  be  pre-implemented  on  a  card. 
FreeMAC  [55]  and  TDMA  MAC  [56]  provide  reconfigurability  to  MAC  pro¬ 
tocol  design.  Their  flexibility  is  however  limited  by  the  coercive  wireless-card 
specifications. 

The  idea  of  using  dedicated  wireless  reconfigurable  platforms  has  also  emerged 
to  address  the  lack  of  reconfigurability  of  existing  hardware.  Software-defined  ra¬ 
dios  (SDRs)  such  as  USRP  [52]  are  a  good  example  of  such  dedicated  platforms. 
USRP  works  along  with  the  GNU  Radio  [57]  framework,  which  offers  flexibility 
thanks  to  its  software-based  processing  and  modular  design  features.  However, 
this  flexibility  is  limited  to  the  PHY  layer.  SORA  [58]  offers  flexibility  along 
with  high  performance  through  a  multi-core  processing  architecture.  However, 
its  design  imposes  difficulties  (in  terms  of  software  complexity)  in  reconfigur¬ 
ing  the  protocol  stack.  Moreover,  platforms  like  AirBlue  [59]  and  WARP  [53] 
transfer  most  of  the  processing  functions  to  Field-Programmable  Gate  Arrays 
(FPGA)  to  improve  performance.  AirBlue  also  has  a  modular  architecture  that 
enables  easier  modifications  in  MAC  design. 

Another  family  of  recent  proposals  are  based  on  the  idea  of  decompos¬ 
ing  wireless  MAC  protocols  into  core  functional  blocks  that  can  be  composed 
through  a  high-level  language  to  enable  reconfigurability  [35,  51,  36,  50].  This 
promising  approach  is  also  at  the  core  of  the  design  of  RcUBe.  In  [35],  the 
architecture  is  built  on  USRP  [57,  52].  The  key  design  idea  is  to  split  MAC 
functionalities  between  the  host  processing  and  the  hardware  (FPGA)  process¬ 
ing  units.  Hence,  time-critical  MAC  functionalities  are  run  on  FPGA  for  better 
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performance  and  timing  control,  while  other  control  functionalities  are  kept  in 
the  host  PC.  Conversely,  in  [51],  all  MAC  primitive  blocks  are  implemented 
on  hardware.  MAC  blocks  are  implemented  (in  WARP  [53])  on  FPGA  boards 
and  a  PowerPC  CPU  is  employed  as  the  main  processing  unit  of  the  archi¬ 
tecture.  In  [36],  the  proposed  architecture  is  built  on  a  resource-constrained 
commodity  WLAN  card.  Unlike  previous  work,  the  authors  propose  the  im¬ 
plementation  of  an  abstract  execution  machine  which,  based  on  conditions  and 
events,  controls  the  execution  of  primitive  building  blocks.  These  blocks  are 
created  by  decomposing  basic  MAC  protocols  and  define  basic  MAC  functions 
(i.e. ,  actions).  Therefore,  MAC  protocols  can  be  flexibly  defined  as  extended 
finite  state  machines  implemented  on  commodity  wireless  cards,  thus  enabling 
real-time  reprogrammability,  without  modifying  the  hardware.  In  a  follow-up 
work  [60] ,  extensions  are  proposed  to  reduce  the  run-time  reconfiguration  delay, 
and  a  control  logic  is  proposed  to  handle  (i.e.,  load,  move,  and  activate)  MAC 
programs.  In  [50],  the  authors  propose  a  programmable  dataplane  architecture 
that  decomposes  wireless  protocols  into  separate  processing  and  decision  planes. 
They  argue  that  their  design  can  work  on  any  multicore  DSP  processor. 

At  the  network  layer,  OpenFlow  [61]  has  been  primarily  proposed  for  testing 
experimental  routing  protocols  on  already  deployed  networks,  particularly  in 
campus  networks.  OpenFlow,  which  is  inspired  by  the  internal  structure  of  an 
Ethernet  Switch,  consists  of  three  main  components,  i.e.,  (i)  an  internal  flow 
table  where  each  flow  is  matched  with  its  corresponding  action,  (ii)  a  secure 
channel  connecting  a  switch  with  an  outer  controller,  and  (iii)  the  OpenFlow 
Protocol  that  offers  a  standard  for  communication  between  an  outer  controller 
and  a  switch.  OpenRoads  or  OpenFlow  Wireless  [62,  63]  are  recent  wireless 
extensions  of  OpenFlow. 

RcUBe  is  similar  in  its  objectives  to  each  of  the  aforementioned  approaches. 
However,  it  departs  from  existing  work  in  several  ways.  First,  RcUBe  enables 
real-time  reconfigurability  at  the  PHY,  MAC,  and  network  layers  of  the  protocol 
stack.  It  is  based  on  the  software-defined  philosophy  at  the  PHY.  The  current 
version  is  based  on  the  GNU  Radio  framework  and  USRP  front-end  to  reach 
a  broader  extent  of  PHY  layer  designs  compared  to  approaches  limited  to  the 
capabilities  of  the  selected  commodity  WLAN  card  [36].  RcUBe  offers  a  strict 
separation  between  decision,  control  and  data  logic  to  obtain  more  flexibility 
compared  to  [50].  As  mentioned  before,  it  uses  a  commercial  and  low-cost  RF 
front-end,  and  also  offers  easier  programmability  compared  to  [51].  Moreover, 
it  offers  real-time  self-optimization  through  its  strictly  separated  design. 

E.  RcUBe  Architecture 

RcUBe  is  structured  into  four  planes,  namely  decision ,  control ,  data,  and  regis¬ 
ter  plane,  as  illustrated  in  Fig.  5.  The  decision  plane  consists  of  user-defined 
decision  algorithms  that  provide  the  control  logic  for  real-time  protocol  opti¬ 
mization,  adaptation,  and  reconfiguration.  In  the  control  plane,  RcUBe  stores 
the  required  logic  for  routing  and  MAC  protocol  execution  as  well  as  data  plane 
management.  Data  processing  takes  place  in  the  data  plane,  while  the  regis- 
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Fig.  5:  RcUBe  Architecture. 


ter  plane  stores  and  manages  access  to  system  parameters  and  environmental 
information.  In  this  section,  each  plane  and  their  interactions  are  discussed  in 
detail. 

F.  Decision  Plane 

The  decision  plane  is  the  component  of  RcUBe  that  provides  decisions,  at  real¬ 
time  if  needed,  for  the  other  planes  based  on  user-defined  decision  algorithms. 
This  plane  is  designed  to  be  optional  upon  design,  meaning  that  it  can  be 
enabled/disabled  for  different  implementations.  The  decision  plane  is  interfaced 
with  the  control  and  data  planes  through  the  register  plane,  which  acts  as  a 
carrier  for  decisions.  This  process  is  described  in  detail  in  Section-I.. 

Decision  Engine.  The  main  component  of  the  decision  plane  is  the  so- 
called  decision  engine ,  which  is  responsible  for  the  implementation  of  user- 
defined  decision  algorithms.  The  decision  engine  contains  a  set  of  user-defined 
algorithms,  each  associated  with  a  sensitivity  list.  Each  item  in  the  sensitiv¬ 
ity  list  describes  conditions  under  which  the  corresponding  decision  algorithm 
is  invoked.  Decision  algorithms  may  vary  from  simple  decision  rules  based  on 
checking  the  value  of  a  variable,  to  more  complex  algorithms  requiring  adaptive 
resource  allocation  as  in  cognitive  radio,  or  even  complex  processing  based  on 
convex  optimization  algorithms. 

The  decision  engine  is  capable  of  executing  multiple  decision  algorithms 
in  parallel.  Algorithms  can  be  executed  in  a  synchronous,  or  asynchronous 
fashion,  as  described  in  detail  in  Section-J..  Decision  algorithms  can  (i)  modify 
a  parameter  in  a  protocol,  (e.g.,  value  of  the  minimum  contention  window  (CW) 
in  a  CSMA  protocol),  (ii)  switch  among  different  modes  within  a  protocol, 
(e.g.,  change  the  modulation  scheme),  (iii)  switch  among  different  protocols 
altogether.  Decisions  are  taken  in  the  decision  plane  by  the  decision  engine, 
but  they  are  not  applied  to  the  system  by  the  decision  engine  itself.  After 
decision  making,  the  decision  engine  writes  the  decision  result  in  the  register 
plane.  Therefore,  decisions  can  be  accessed  by  the  data  and  control  planes  to  be 
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adopted  in  their  execution  logics.  This  will  be  further  explained  in  the  following 
sections. 

The  execution  logic  of  the  decision  engine  can  be  illustrated  with  the  fol¬ 
lowing  toy  examples.  First,  consider  a  decision  algorithm  that  is  defined  to 
decide  the  modulation  scheme  of  the  PHY  layer  according  to  the  current  SINR 
value.  In  this  case,  SINR  values,  obtained  by  the  data  plane,  and  user-defined 
threshold  value(s),  which  are  stored  in  the  register  plane,  are  included  in  the 
sensitivity  list  and  monitored  by  the  decision  plane.  When  the  measured  SINR 
value  decreases  below  the  level  of  the  user-defined  threshold  value,  the  decision 
engine  triggers  the  execution  of  the  decision  algorithm  to  select  the  appropriate 
modulation  scheme.  The  selected  modulation  scheme  parameter  is  afterwards 
written  into  the  register  plane  and  thus  becomes  accessible  by  the  data  plane. 
Subsequently,  the  modulation  scheme  parameter  will  be  accessed  by  the  data 
plane  to  adapt  execution  at  the  physical  layer  accordingly.  Second,  consider  a 
network  with  centralized  control.  The  central  node  (access  point  or  base  sta¬ 
tion)  runs  a  decision  algorithm  that  adaptively  selects  a  specific  MAC  protocol 
and  imposes  it  to  the  client  nodes  through  control  packets,  based  on  the  offered 
load  to  the  network.  For  example,  when  the  offered  load  is  above  a  certain 
value,  Time-Division  Multiple  Access  (TDMA)  based  protocols  like  TRAMA 
[64]  are  selected,  since  they  give  better  results  in  energy  efficiency  and  channel 
utilization  when  compared  to  Carrier  Sense  Multiple  Access  (CSMA)  based  al¬ 
gorithms  [64].  The  process  is  executed  as  in  the  previous  example,  and  at  the 
last  step,  the  decision  is  written  into  a  parameter  that  indicates  MAC  protocol 
change.  The  decision  controller  of  the  control  plane,  which  will  be  explained  in 
the  following  section,  detects  the  parameter  change  and  loads  the  newly  selected 
MAC  protocol  design  from  the  MAC  look- up-table  (LUT)  library  of  the  register 
plane  to  the  MAC  execution  engine  of  the  control  plane. 

G.  Control  Plane 

The  control  plane  is  designed  to  handle  the  control  logic  of  data  processing, 
which  takes  place  in  the  data  plane.  The  control  plane  decides  the  sequence  of 
data  processing  functionalities  that  will  be  executed  in  the  data  plane.  It  does 
this  mainly  by  incorporating  two  execution  engines.  The  two  execution  engines 
contain  and  execute  a  FSM  implementing  the  logic  of  either  the  MAC  or  the 
network-layer  protocol  execution,  respectively.  The  control  plane  further  con¬ 
tains  a  design  controller  to  manage  the  logic  executed  in  the  execution  engines, 
as  shown  in  Fig. 5. 

Execution  Engines.  Execution  engines  are  designed  to  execute  FSMs  that 
define  MAC  and  routing  protocols.  The  FSMs  are  designed  as  Extended  Finite 
State  Machines.  Each  FSM  contains  symbolic-states  and  an  extended  state 
transition  definition.  Each  state  transition  is  designed  as  a  triplet  that  contains 
events ,  conditions ,  and  actions.  Each  state  transition  occurs  if  events  take  place 
while  conditions  are  fulfilled.  As  a  result  of  that,  the  corresponding  actions  are 
invoked.  Actions  are  mainly  designed  to  invoke  a  data  processing  functionality 
in  the  data  plane,  while  events  may  occur  as  a  result  of  data  processing  actions. 
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(b)  Transmitter  path  of  the  FSM. 


Fig.  6:  Sample  FSM  architecture  taking  place  in  the  MAC  execution  engine  at 
the  control  plane.  State  transitions  are  based  on  the  triplet  (events,  conditions, 
actions) . 


Conditions  are  controlled  via  the  register  plane,  where  system  parameters  and 
environmental  variables  are  stored.  To  simplify  the  definition  of  new  protocols, 
the  state  transitions  of  FSMs  are  defined  by  LUTs,  which  are  stored  in  the  reg¬ 
ister  plane  to  create  a  library  of  different  protocols,  and  loaded  on  the  execution 
engines  to  be  transformed  into  FSMs  for  execution.  A  sample  list  of  events,  con¬ 
ditions,  and  actions  are  shown  in  Table- 1.  Each  component  and  its  properties 
are  explained  in  detail  in  the  following  sections.  The  core  design  idea  of  using 
FSM  structures  is  adopted  from  [36].  Here,  it  is  reformulated  at  a  different  level 
of  abstraction  to  fit  our  design  and  to  include  both  MAC  and  routing  protocols. 
Figure  6  depicts  a  sample  FSM  of  a  MAC  layer  design  described  for  the  CR 
scenario  at  [43]. 

Design  Controller.  This  component  is  responsible  for  control  and  man¬ 
agement  of  the  process  of  loading  routing  and  MAC  layer  protocols,  which  are 
defined  by  the  user  as  LUTs  in  the  register  plane,  into  the  execution  engines.  It 
also  controls  switching  between  different  protocols  (pre-defined  in  the  libraries 
at  the  register  plane)  on-the-fly  based  on  the  decisions  taken  by  the  decision 
plane.  The  design  controller  monitors  parameters  in  the  register  plane  that 
identify  which  protocols  need  to  be  executed  by  the  execution  engines.  These 
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parameters  can  either  be  updated  in  real  time  by  decisions  taken  in  decision 
plane,  or  offline  by  the  system  designer. 


Table  1:  List  of  Events,  Conditions,  and  Actions. 


Event 

Condition 

Action 

mac_has_data 

channel !  =busy 

set  Timer  () 

channel-busy 

backoff  !=0 

resume_bk() 

timer_end 

backoff==0 

set_bk() 

bk_end 

RTS_lim!=lim 

freeze_bk() 

bk_end 

RTS_lim==lim 

send_RTS() 

rts_sent 

num_of_rtran! =lim 

rec_RTS() 

dts_sent 

mini  of  rtran==lim 

send_CTS() 

cts_sent 

sch.queue ! =empty 

rec_CTS() 

ack_sent 

sch_queue==empty 

send_DTS() 

channeLidle 

queues  !=  empty 

set  _opt  imizat  ion  .enable  ( ) 

data_sent 

rec_DTS() 

cts_rec 

send_data() 

cts_not  jrec 

rec.dataQ 

ack_rec 

send_ACK() 

ack_not_rec 

rec_ACK() 

cts_not_rec 

send_next_data() 

ack_not_rec 

send_data() 

ack_not_rec 

recdRTSQ 

data_not_rec 

send_next_data() 

mac_data_req 

send_data_mac() 

optimization_end 

set_rout_has_data() 

mac_has_data 

rec_data_mac() 

mac datajrec 

identify data() 

H.  Data  Plane 

The  data  plane  is  where  data  processing  and  data/control  queues  are  handled.  It 
also  contains  two  libraries  of  pre-defined  primitive  building  blocks,  actions  and 
events ,  and  interacts  with  the  RF  front-end/USRP.  Briefly,  data  coming  from 
upper  layers  are  received  at  the  network  layer,  and  after  header  formation  at 
the  network  and  MAC  layer,  are  transformed  into  digital  waveforms  at  the  PHY 
layer,  which  is  implemented  in  GNU  Radio.  Subsequently,  digital  waveforms  are 
converted  into  analog  waveforms  by  the  digital-to-analog  converter  (DAC)  in  the 
RF  front-end. 

The  data  plane  has  three  main  design  features  that  separate  it  from  tra¬ 
ditional  systems.  First,  it  decomposes  protocols  into  primitive  building  blocks. 
In  this  way,  the  data  plane  has  a  modular  and  flexible  structure  and  can  be 
reconfigured  by  changing  the  sequence  of  execution  of  these  primitive  building 
blocks.  Primitive  building  blocks  include  actions  and  events.  Actions  are  func- 
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tional  blocks  that  execute  data  processing,  whereas  events  are  indicator  blocks 
for  acquiring  results  of  data  processing.  Second,  as  briefly  discussed  in  Section- 
G.,  is  the  separation  between  data  processing  and  its  control  logic.  This  means 
that  the  control  plane  only  decides  the  sequence  of  data  processing  functional¬ 
ities  with  its  executing  engines,  while  data  processing  itself  takes  place  in  the 
data  plane.  Thanks  to  this  separation,  RcUBe  is  able  to  define  different  pro¬ 
tocols  just  by  changing  the  sequence  of  data  processing  actions.  For  instance, 
suppose  the  control  plane  decides  to  send  an  RTS  packet  to  a  neighbor  node  by 
invoking  the  action  named  “send_RTS()”.  After  this,  the  data  plane  acquires 
the  next  hop  information  from  the  routing  table  stored  in  the  register  plane  and 
accordingly  forms  the  RTS  packet  structure  at  the  MAC  layer  and  transmits 
it  at  PHY  layer  by  using  the  corresponding  GNU  Radio  flowgraph.  Note  that 
all  the  parameters  of  protocols  at  different  layers  are  stored  as  variables  in  the 
register  plane.  Therefore,  any  parameter  of  any  protocol  at  any  layer  can  be 
reconfigured  on-thc-fly  through  the  register  plane.  Protocol  variables  can  be 
accessed  “intrinsically”  by  the  decision  plane  as  a  result  of  a  decision  algorithm, 
or  by  the  data  plane  according  to  specific  data  acquisitions,  or  “extrinsically” 
by  the  protocol  designer. 

Actions.  Actions  are  primitive  building  blocks  formed  as  a  decomposition  of 
MAC  and  network  layer  protocols.  Actions  are  mainly  designed  to  handle  data 
plane  processes.  Nevertheless,  they  are  also  used  to  update  parameters  stored 
in  the  register  plane  and  to  provide  synchronization  between  MAC  and  routing 
FSMs,  particularly  to  exchange  data  between  each  other.  They  can  be  classified 
based  on  their  functionalities  as  (i)  actions  to  transmit /receive  data  and  control 
packets  (e.g.,  ACK,  CTS,  RTS,  among  others),  including  processing  of  packets 
with  headers  and  selecting  corresponding  PHY  flowgraph  as  subsystems,  (ii) 
actions  to  reconfigure  parameters  of  protocols  at  different  layer,  (iii)  timing 
functionalities,  (iv)  routing  functionalities  (e.g.,  forwarding,  exchanging  route 
discovery  messages),  (v)  data  transfer  between  different  layers. 

Events.  Events  are  primitive  building  blocks  that  are  responsible  for  indi¬ 
cating  the  results  of  data  processing  functions,  so  they  maintain  state  transitions 
and  protocol  executions,  as  discussed  in  Section-G.. 

To  handle  MAC  and  routing  layer  processing,  we  have  built  our  primitive 
building  blocks  as  actions.  Moreover,  for  routing  protocols,  we  have  defined  a 
routing  table  in  the  register  plane.  For  PHY  layer  processing,  since  our  design 
was  developed  on  USRP,  we  have  taken  advantage  of  the  open  source  GNU  Ra¬ 
dio  framework.  GNU  Radio  provides  an  extensive  library  of  signal  processing 
blocks.  Although  many  blocks  already  exist,  new  blocks  can  be  added  if  needed. 
Thus,  one  can  define  a  PHY  protocol  simply  by  connecting  these  blocks  so  as 
to  build  a  flowgraph.  Thanks  to  this  granularity,  our  framework  can  enable 
reconfigurability  in  three  main  ways.  First,  according  to  the  PHY-layer  prefer¬ 
ences  of  the  MAC  protocol,  the  control  plane  can  switch  between  pre-defined 
flowgraphs.  Second,  the  parameters  of  blocks  that  form  flowgraphs  can  be  up¬ 
dated  through  the  register  plane  based  on  the  decision  engine.  Finally,  users 
can  simply  connect  different  blocks  to  implement  new  PHY  protocols. 
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I.  Register  Plane 

The  register  plane  is  responsible  for  saving  both  environmental  variables/readings 
(e.g.,  SINR)  and  state  variables  at  the  different  layers  (e.g.,  modulation,  trans¬ 
mit  power,  size  of  minimum  CW,  among  others).  It  acts  as  a  reference  plane 
for  the  other  planes  (i.e.,  decision,  data,  control),  and  it  can  be  updated  by 
both  the  decision  plane  (based  on  results  from  executed  algorithms)  and  by 
the  data  plane  according  to  environmental  sensed  data  or  incoming  information 
from  neighboring  nodes.  It  can  be  accessed  by  all  planes.  The  register  plane  also 
hosts  two  libraries  of  LUTs  that  encapsulate  MAC  and  routing  protocols.  These 
libraries  are  designed  to  store  custom  MAC  and  routing  protocols.  Therefore, 
users  can  easily  create  their  own  cross-layer  protocol  libraries,  thus  achieving 
fast  prototyping.  Moreover,  the  control  plane  may  use  these  libraries  to  switch 
protocol  on-the-fly  based  on  decisions  taken  at  the  decision  plane.  The  regis¬ 
ter  plane  also  hosts  neighbor  and  routing  tables,  which  may  be  accessed  and 
updated  by  the  decision  and  data  planes  depending  on  specific  protocol  designs. 

J.  Architectural  Considerations 

Delay.  The  current  version  of  the  RcUBe  framework  is  fully  implemented 
in  software  and  it  mainly  runs  on  the  host  PC.  This  design  choice  was  made 
to  fully  leverage  the  flexibility  of  the  GNU  Radio  framework  for  PHY  layer 
processing.  As  discussed  and  observed  before,  pure  software  implementations 
of  the  protocol  stack  offer  great  flexibility  at  the  expense  of  higher  processing 
latency  and  coarser  control  on  the  timing  of  operations  [35,  65].  In  Section 
L.,  we  also  observed  delays  that  prevent  our  framework  from  meeting  time- 
critical  deadlines  of  the  implemented  MAC  protocols.  Thus,  as  a  temporary 
fix,  we  defined  and  adopted  specialized  guard  times  to  maintain  the  timing 
and  execution  of  the  MAC  protocols.  The  software-hardware  implementation 
boundaries  of  our  framework  are  flexible,  and  we  are  currently  working  to  shift 
portions  of  the  execution  engines  to  a  soft-core  processor  implemented  on  FPGA, 
thus  combining  the  flexibility  offered  by  high-level  programming  languages  with 
the  processing  efficiency  of  hardware  implementations. 

Control  and  Data  Packet  Handlers.  RcUBe  offers  data  and  control 
packet  handlers  through  actions  defined  in  the  control  plane.  The  functionalities 
of  these  actions  are  two-fold.  Firstly,  they  provide  switching  between  different 
physical  layer  receiver  flowgraphs  depending  on  the  packet  type  (i.e.,  control  or 
data  packet).  Secondly,  they  are  responsible  for  the  implementation  of  MAC 
layer  functionalities.  For  example,  consider  the  action  named  “rec-RTS()”  . 
This  action  is  defined  for  handling  both  in  physical  and  MAC  layers,  control 
packets  named  Request-to-Send  (RTS). 

Asynchronous  or  Synchronous  Decisions.  RcUBe  supports  both  asyn¬ 
chronous  and  synchronous  decisions.  A  decision  algorithm  can  be  either  ex¬ 
ecuted  synchronously,  i.e.,  triggered  by  another  plane,  or  asynchronously.  In 
particular,  for  synchronous  execution,  a  flag  variable  should  be  defined  into  the 
sensitivity  list  of  the  corresponding  algorithm  in  the  decision  plane.  Depending 
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on  which  protocol  is  adopted  in  the  register  plane  the  activation  time  of  this  flag 
by  the  control  plane  is  affected.  After  execution  of  the  corresponding  decision 
algorithm,  another  flag  variable  should  be  enabled  to  inform  the  decision  plane 
about  completion  of  execution.  However,  sometimes  the  synchronous  approach 
may  be  harmful  for  the  protocol  stack  performance  because  of  the  delays  asso¬ 
ciated  with  execution  of  decision  algorithm.  In  these  cases,  decision  algorithms 
can  be  executed  asynchronously  without  expecting  an  acknowledgement  from 
the  decision  plane.  In  this  way,  the  rest  of  the  system  does  not  need  to  wait  for 
the  optimized  parameters  obtained  by  the  decision  algorithm  to  execute  the  pro¬ 
tocol  stack.  Therefore,  one  can  either  choose  to  use  asynchronous  decisions  to 
meet  hard  deadlines  in  protocol  stack  execution  and  operate  with  non-optimized 
parameters,  or  wait  for  optimized  parameters  at  the  price  of  extra  delay. 

Centralized  and  Distributed  Control.  RcUBe  can  be  used  in  networks 
with  both  centralized  and  distributed  control.  In  Section-K.,  we  will  discuss  the 
implementation  of  a  complex  protocol  in  a  network  with  distributed  control.  For 
networks  with  centralized  control,  the  following  toy  example  can  be  explanatory. 
Consider  a  cluster-based  network  with  TDMA  scheduling.  Here,  all  decisions, 
e.g.,  number  of  time  slots  to  be  allocated  to  each  user,  are  taken  by  the  cluster 
head  and  forwarded  to  cluster  members  through  control  packets.  At  cluster 
members,  RcUBe  can  be  configured  with  a  disabled  decision  plane.  In  this 
way,  children  nodes  can  only  decode  received  control  packets  and  write  enclosed 
information  (number  of  slots  assigned  to  them  for  next  session)  to  their  register 
plane.  At  the  cluster  head,  an  algorithm  running  in  the  decision  plane  assigns 
time  slots  to  different  nodes. 

K.  Example  Design  Implementation 

Our  primary  objective  is  to  demonstrate  the  flexibility,  real-time  reconfigura¬ 
bility,  and  self-optimization  capabilities  offered  by  RcUBe.  To  this  end,  we 
implemented  within  the  RcUBe  framework  a  joint  ROuting  and  Spectrum  Al¬ 
location  algorithm  (ROSA)  algorithm  [43].  The  ROSA  algorithm  was  chosen 
because  of  its  ability  to  demonstrate  and  take  advantage  of  all  the  proposed 
features  of  the  RcUBe  framework.  ROSA  is  based  on  a  cross-layer  architecture 
(PHY,  MAC,  and  Routing)  with  complex  interactions  and  employs  a  decen¬ 
tralized  optimization  algorithm.  Moreover,  it  continuously  relies  on  real-time 
decisions  and  reconfigurations  as  it  adapts  its  PHY,  MAC,  and  routing  behavior 
according  to  time-varying  traffic  demands,  network  topology,  and  interference 
profile.  Additionally,  on  the  control  channel  ROSA  employs  an  IEEE  802.11- 
like  MAC  protocol,  thus  giving  us  an  opportunity  to  prove  that  RcUBe  can 
support  both  conventional,  current  state-of-the-art  protocol  designs  as  well  as 
newly  proposed,  novel  cross-layer  protocols. 

Definition  of  the  Algorithm.  ROSA  [43]  is  a  CR  algorithm  that  tries  to 
maximize  the  network  throughput  through  cross-layer  jointly  optimal  routing, 
dynamic  spectrum  allocation,  scheduling,  transmit  power  control  and  interfer¬ 
ence  avoiding  waveform  selection. 

The  algorithm  was  developed  specifically  for  secondary  users  of  a  cognitive 
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Fig.  7:  Routing  execution  engine  for  ROSA. 


ad  hoc  networks,  which  send  their  data  opportunistically  without  deteriorating 
on-going  (primary)  communications.  ROSA  is  based  on  two  different  chan¬ 
nels.  Secondary  users  contend  for  spectrum  access  and  exchange  local  control 
information  on  the  common  control  channel  (CCC),  while  the  data  channel 
(DC)  consists  of  equally-sized  discrete  frequency  bands  to  be  allocated  by  users. 
ROSA  employs  an  IEEE  802.11-like  protocol  with  three-way  handshaking  and 
spectrum  sensing  for  enabling  Collaborative  Virtual  Sensing  (CVS)  [43].  Each 
node  can  acquire  local  information  through  explicit  exchange  of  control  mes¬ 
sages,  by  overhearing  control  packets,  or  through  spectrum  sensing.  Therefore, 
each  backlogged  node  is  able  to  acquire  spectrum  and  queuing  information  from 
its  neighbors  and  can  accordingly  select  the  next  hop  and  spectrum  band  to 
be  used  for  transmission.  An  additional  control  packet,  named  Data  Transmis¬ 
sion  reservation  (DTS),  was  defined,  in  addition  to  Request-to-Send  (RTS)  and 
Clear-to-Send  (CTS),  to  be  responsible  for  informing  the  neighbor  nodes  about 
spectrum  reservation  and  transmit  power. 


(a)  Testbed  deployment  in  EE 
Dept. 


(b)  Schematic  overview  of  the  5-node 
testbed  setup. 


Fig.  8:  A  view  from  our  testbed  implementation. 

ROSA  is  based  on  a  backpressure  principle  and  performs  a  joint  routing, 
dynamic  spectrum  allocation,  and  waveform  selection  algorithm  that  is  triggered 
whenever  a  backlogged  node  senses  the  CCC  to  be  idle.  Each  node,  source 
or  intermediate  node,  maintains  a  separate  queue  for  each  session,  defined  by 
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a  fixed  node-destination  pair.  The  algorithm  identifies  a  set  of  feasible  next 
hops  for  each  session  and  for  each  of  them  solves  a  link  capacity  maximization 
problem  to  determine  optimal  spectrum  and  power  allocations.  Subsequently, 
it  calculates  a  so-called  spectrum  utility  function,  defined  as  the  product  of 
maximum  link  capacity  and  maximum  differential  backlog  on  that  link.  It  then 
selects  session  and  link  (next  hop)  with  maximum  spectrum  utility  function. 
According  to  the  link  chosen,  it  decides  on  a  modulation  scheme  (i.e. ,  BPSK, 
QPSK)  based  on  the  current  estimated  SINR  on  the  chosen  channel  at  the 
transmitter  node.  Finally,  the  node  uses  the  relative  spectrum  utility  function 
to  select  the  size  of  the  CW,  which  determines  the  priority  of  being  scheduled 
for  transmission  as  compared  to  neighboring  competing  nodes. 

Interpretation  of  the  Algorithm.  The  ROSA  algorithm  and  its  func¬ 
tionalities  are  implemented  in  the  RcUBe  architecture  as  follows.  The  control 
plane  contains  the  execution  engines  for  routing  and  MAC  protocols  that  are 
illustrated  in  Fig.  7  and  Fig.  6  respectively.  The  routing  protocol  FSM  is  de¬ 
signed  in  such  a  way  as  to  trigger  the  decision  algorithm  that  takes  place  in  the 
decision  plane.  Moreover,  it  also  invokes  actions  to  form  and  handle  network 
layer  headers,  in  coordination  with  the  MAC-layer  FSM,  through  the  routing 
information  contained  in  the  routing  table  in  the  register  plane.  Therefore,  the 
routing  FSM  prepares  packets  to  be  sent  to  the  MAC  layer  or  identifies  received 
packets  from  the  MAC  layer.  The  FSM  design  of  the  MAC  execution  engine 
is  shown  in  Fig.  6  with  two  subfigures  depicting  receiver  and  transmitter  path 
that  are  connected  to  each  other  through  the  shared  state  IDLE.  This  FSM 
is  formed  based  on  common  primitive  building  blocks  that  can  be  used  for  any 
802.11-like  MAC  protocol. 

The  decision  plane  hosts  and  executes  the  decision  algorithm  of  ROSA  that 
decides  next  hop,  spectrum  portion  to  be  occupied,  waveform  and  size  of  con¬ 
tention  window.  In  order  to  trigger  the  execution  of  this  algorithm,  the  decision 
engine  monitors  an  enable  parameter  (  “optimization-enable”)  that  is  defined  in 
its  sensitivity  list.  This  flag  is  raised  as  a  consequence  of  an  action  invoked  by 
the  control  plane  (routing  execution  engine)  and  stored  in  the  register  plane.  Af¬ 
ter  execution  of  the  decision  algorithm,  the  updated,  optimized  parameters  are 
written  in  the  register  plane.  Moreover,  an  event  is  raised  so  as  to  indicate  the 
end  of  the  decision  algorithm  execution.  This  event  is  used  as  a  state-transition 
trigger  in  the  routing  execution  engine.  This  process  is  a  good  example  of  syn¬ 
chronous  decisions,  as  discussed  in  Section  J..  The  data  plane  contains  the  data 
processing  of  PHY,  routing  layers  and  queues  (sessions)  for  data  packets.  The 
register  plane  is  the  storage  space  for  all  the  system  and  environmental  param¬ 
eters  that  can  be  accessed  by  other  planes.  Apart  from  that,  it  also  hosts  the 
neighbor  table  that  is  pre-defined  for  the  ROSA  algorithm,  routing  table,  and 
LUT  libraries  for  other  MAC  and  routing  protocols. 

L.  Experimental  Evaluation 

In  this  section,  we  present  experimental  results  obtained  by  implementing  the 
ROSA  protocol  suite  in  RcUBe  in  our  testbed  depicted  in  Fig.  8-a.  Each  wireless 
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(a)  Carrier  Frequency  for  Dynamic  and  Fixed  Allocations. 


SINR  for  Dynamic  and  Fixed  Allocations 


(b)  SINR  for  Dynamic  and  Fixed  Allocations. 


(c)  Session  Backlog  for  Dynamic  and  Fixed  Allocations. 


Fig.  9:  Testbed  Results  for  2-Node  Setup. 


software-defined-radio  (SDR)  node  consists  of  a  Linux-PC,  hosting  the  RcUBe 
framework,  connected  via  Gigabit  Ethernet  (GigE)  with  a  USRPN  —  210  SDR 
device,  developed  by  Ettus  [52] .  Each  PC  is  equipped  with  an  Intel  i7  —  3610QM 
2.30  GHz  processor,  8  GByte  RAM.  The  operating  system  is  Ubuntu  Precise 
Pangolin  12.04.  We  used  up  to  5  nodes  named  Nl,  N2,  N3,  N4,  N5,  while 

39 


Approved  for  Public  Release;  Distribution  Unlimited 


an  extra  6th  node  (named  TR)  was  used  for  capturing  trace  data  for  offline 
processing. 

We  used  an  orthogonal  frequency-division  multiplexing  (OFDM)  scheme 
with  200  subcarriers  at  the  PHY  layer.  Each  OFDM  subcarrier  is  modulated  ei¬ 
ther  with  Binary-Phase-Shift-Keying  (BPSK)  or  Quadrature-Phase-Shift-Keying 
(QPSK).  In  each  communication  link,  the  modulation  scheme  is  picked  accord¬ 
ing  to  a  decision  made  by  the  algorithm  based  on  the  SINR  estimation  of  this 
link.  Specifically,  we  defined  a  SINR  threshold  at  9  dB.  Therefore,  if  the  SINR 
estimate  of  the  chosen  link  is  lower  than  this  threshold,  the  algorithm  selects 
BPSK,  otherwise  it  selects  QPSK.  We  considered  a  CCC,  which  has  a  center 
frequency  at  2.41  GHz  with  a  bandwidth  of  1  MHz.  We  considered  a  DC  with 
five  channels  with  center  frequencies  at  2.418,  2.426,  2.434,  2.442,  and  2.450GHz 
respectively,  and  1  MHz  bandwidth. 

We  considered  two  networks  of  different  size,  i.e.,  2-Node  and  5-Node  setups, 
to  focus  on  different  features  and  capabilities  of  the  framework. 

2-Node  Setup.  The  first  set  of  experiments  were  implemented  on  a  2-Node 
setup,  where  there  was  1-session  defined  from  N1  to  N2.  In  these  experiments, 
our  intention  was  to  “degrade”  the  ROSA  decision  algorithm  to  a  dynamic  spec¬ 
trum  allocation  and  optimal  waveform  selection  algorithm  only,  to  concentrate 
on  PHY  reconfigurations.  Since  we  had  only  one  hop,  the  differential  backlog 
factor  was  not  considered  and  the  ROSA  algorithm  simply  determined  the  best 
spectrum  allocation  and  waveform  selection  to  maximize  the  channel  capacity. 
In  this  set  of  experiments,  we  considered  two  different  scenarios.  In  the  first  one 
we  used  dynamic  spectrum  allocation  as  in  ROSA,  while  in  the  second  we  used 
fixed  spectrum  allocation.  The  algorithm  just  decides  the  modulation  scheme 
to  be  picked  based  on  the  SINR  estimate.  Each  experiment  lasted  for  60  s. 

Figure  9  illustrates  the  experimental  results  acquired  by  the  setup  described 
above.  Figure  9- (a)  depicts  the  frequency  allocation  results  for  both  cases  (dy¬ 
namic  and  fixed  allocation),  while  Fig.  9-(b)  illustrates  the  corresponding  SINR 
estimates.  It  is  clear  that  dynamic  spectrum  allocation  obtains  a  higher  SINR 
profile  when  compared  to  fixed  spectrum  allocation.  As  a  result,  dynamic  spec¬ 
trum  allocation  can  obtain  higher  transmission  rates  as  it  can  select  higher-order 
modulation  schemes  (i.e.,  QPSK  instead  of  BPSK).  This  effect  can  be  observed 
in  Fig  9-(c),  which  shows  the  size  of  session  backlog  for  both  cases.  It  can  be 
concluded  that  when  spectrum  is  allocated  dynamically,  session  backlog  gets 
cleared  faster  than  with  fixed  spectrum  allocation.  Overall,  real-time  decisions 
about  the  PHY  layer  parameters  (i.e.,  modulation  and  carrier  frequency  selec¬ 
tion)  can  be  applied  in  real  time. 

5-Node  Setup  The  second  set  of  experiments  was  implemented  on  a  5- 
Node  setup,  where  there  was  1-session  defined  from  N1-N5,  active  for  120  s. 
Based  on  ROSA,  each  node  has  a  pre-defined  neighborhood  table  for  routing 
purposes.  We  defined  the  neighborhood  in  such  a  way  that  N2,  N3  and  N4  were 
intermediate  nodes,  only  one  hop  away  from  destination  node  N5.  Each  node 
can  overhear  control  packets  from  neighbors  through  the  CCC.  The  testbed 
schematic  in  Fig.  8-(b)  shows  the  neighborhood  connections  between  the  nodes. 

Figures  10  and  11  depict  the  results  acquired  from  the  experiments  conducted 
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(a)  Differential  Backlog  for  Links  [N1-N2],  [N1-N3],  and  [Nl- 
N4], 


Utility  function  for  Links  [N1-N2],  [N1-N3J,  and  [N1-N4] 


Backlogged  Queue  Sizes  of  N2,  N3,  N4 


Fig.  10:  Testbed  Results  for  5-Node  Setup  with  1-Session  at  [N1-N5]. 


in  the  5-Node  setup.  Figure  10-(a)  illustrates  the  differential  backlog  for  the  links 
between  the  source  and  three  intermediate  nodes.  During  the  first  10  seconds, 
N2  gets  most  of  the  traffic  load  from  N1  as  its  link  has  higher  SINR  compared  to 
N3  and  N4.  This  can  also  be  observed  in  Fig.  11.  Accordingly,  higher  spectrum 
utility  is  achieved  because  of  the  higher  channel  capacity.  Moreover,  since  link 
N1-N2  has  high  SINR  values,  it  uses  QPSK  modulation,  which  results  in  higher 
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(a)  Best  SINR  values  for  links  [N2-N5],  [N3-N5],  and  [N4-N5] 
after  dynamic  spectrum  selection. 


Best  SINR  values  for  links  [N1-N2],  [N1-N3],  and  [N1-N4]  after  dynamic  spectrum  allocation 


(b)  Best  SINR  values  for  links  [N1-N2],  [N1-N3],  and  [N1-N4] 
after  dynamic  spectrum  selection. 


Fig.  11:  Best  SINR  values  for  links  in  the  5-Node  Setup  with  1-Session  at  [Nl- 
N5]. 


transmission  rates.  Between  the  10th  and  20th  second,  consequent  to  its  fast 
start,  N1-N2  reaches  a  lower  differential  backlog  value  compared  to  the  other 
links,  because  of  backlogged  packets  at  N2.  Therefore,  the  utility  functions  of 
the  two  links  are  balanced,  as  can  be  observed  from  Fig.  10- (b),  and  eventually 
the  traffic  is  shared  equally  between  the  intermediate  nodes.  After  the  20th 
second,  it  can  be  observed  that  each  one  hop  link  has  similar  utility  function, 
which  makes  them  receive  incoming  traffic  from  the  source  equally.  However, 
it  can  be  observed  from  Fig.  10-(a),  that  the  differential  backlog  value  of  the 
N1-N2  link  reaches  the  level  of  others  around  the  120th  second.  The  reason 
behind  this  can  be  understood  by  observing  Fig.  10-(c),  where  the  backlog 
queue  sizes  of  intermediate  nodes  are  shown.  The  high  traffic  offered  from  the 
N1-N2  link  and  the  low  SINR  values  at  the  N2-N5  link  cause  N2  to  become 
highly  backlogged.  However,  after  the  20th  second,  as  N2  is  more  likely  to  be 
scheduled  for  transmission  due  to  the  narrow  contention  window  selection  (which 
is  directly  related  to  the  utility  function  of  the  link)  the  number  of  backlogged 
packets  decreases.  Therefore,  the  N1-N2  link  has  a  higher  differential  backlog 
than  before,  even  if  it  did  not  get  more  traffic. 

It  can  be  concluded  that  each  node  was  capable  of  taking  decisions  and 
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reconfigure  in  real-time  PHY  (i.e. ,  waveform  selection,  spectrum  allocation), 
MAC  (i.e.,  size  of  CW)  and  network  layer  parameters  (i.e.,  next  hop  decisions) 
through  a  decision  algorithm. 
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M.  Software  Defined  Radio  All-spectrum  Cognitive  Chan¬ 
nelization 

Cognitive  Radio  (CR)  has  emerged  as  a  promising  technology  for  improving 
spectrum  utilization  efficiency  [66,  67,  68,  69].  CR  systems  have  been  envisioned 
as  intelligent,  flexible,  radio  systems  capable  of  reconfiguring  their  own  parame¬ 
ters  (i.e.,  modulation,  waveform,  carrier  frequency  etc.)  in  an  autonomous  way 
[70],  and  may  compare  favorably  to  traditional  communication  systems  which 
are  missing  such  capabilities  and  often  criticized  for  wasting  resources.  The 
evolution  of  software  defined  radio  (SDR)  technology  enables  implementation 
of  such  “smart”  radios,  capable  of  reconfiguring  crucial  communication  param¬ 
eters  in  software,  offering  that  way  extreme  agility  over  hardware  radios  with 
fixed  characteristics.  CR  users,  also  called  secondary  users,  may  coexist  with 
primary  users  (spectrum  licensees)  and  be  allowed  to  access  licensed  spectrum 
provided  that  their  interference  to  primary  users  is  maintained  below  a  prespec¬ 
ified  threshold  (transmission  in  grey  spaces  [72,  71,  73,  75,  74]).  Alternatively, 
CR  users  may  transmit  in  holes  (white  spaces)  of  unutilized  spectrum  which 
may  not  be  the  preferred  accessing  method  in  highly  occupied  spectrum  bands 
[76].  SDR  implementation  of  cognitive  radios  has  been  considered  in  [78,  79], 
where  secondary  users  sense  the  spectrum  and  transmit  in  spectrum  holes  that 
may  or  may  not  be  contiguous. 

Over  the  past  several  years  we  have  been  working  on  and  contributing  to 
cognitive  radio  technologies,  where  primary  and  secondary  users  coexist  and 
utilize  the  same  spectrum  concurrently.  In  this  project,  we  consider  secondary 
links  coexisting  with  unknown  primary  users  and  we  experimentally  evaluate 
the  CR  spectrum  access  methodology  proposed  in  [72,  71,  73,  75,  74]  on  a 
low-cost,  indoors  deployed  SDR  testbed.  Spectrum  access  is  achieved  through 
adaptive  channelization  where  adaptively  optimized  waveforms  are  assigned  to 
secondary  users.  Waveforms  are  designed  to  maximize  the  output  signal-to- 
interference  plus  noise  ratio  (SINR)  at  the  output  of  the  secondary  receivers, 
and  at  the  same  time,  maintain  the  interference  to  the  primary  users  below  a 
prespecified  threshold  [72,  71,  73,  75,  74,  77]. 

Our  testbed  is  deployed  in  an  RF-cluttering  lab  environment,  assumes  un¬ 
known  primary  users  and  utilizes  low-cost,  commodity  transceivers  whose  inter¬ 
nal  oscillators  are  quite  unstable,  and  thus  introduce  carrier  frequency  offsets 
(CFOs)  between  any  transmitter-receiver  pair.  Thus,  the  waveform  design  prob¬ 
lem  presented  in  [73] ,  has  to  be  coupled  with  estimation  of  the  unknown  multi- 
path  channel  coefficients  and  carrier  frequency  offsets.  In  other  words,  while  our 
main  objective  in  this  project  is  the  experimental  evaluation  of  the  CR  access 
methodology  in  [72,  71,  73,  75,  74],  as  a  side  theoretical  result  we  develop  a  tech¬ 
nique  for  blind  estimation  of  multipath  channel  coefficients  and  CFOs.  More 
specifically,  we  propose  a  generalized  version  of  the  subspace-based  estimation 
procedure  presented  in  [80].  Then,  using  three  commercial  SDR  transceivers,  we 
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experimentally  validate  the  concept  of  cognitive  channelization  and  spectrum 
sharing  between  primary  and  secondary  users,  in  terms  of  instantaneous  SINR 
and  bit-error-rate  (BER). 

N.  System  Model 

In  this  work,  we  consider  a  low-cost  secondary  link  (secondary  transmitter- 
receiver  pair)  operating  in  the  same  frequency  as  a  primary  user.  Both  the 
secondary  and  primary  users  are  assumed  to  transmit  binary  antipodal  infor¬ 
mation  symbols  bk(i)  £  {— 1, 1},  i  =  0, . . . ,  J  —  1  and  k  =  1,  2  at  a  rate  1/T, 
modulated  by  a  signal  waveform  dk(t)  of  duration  T.  If  fc  denotes  the  com¬ 
mon  carrier  frequency,  then  the  transmitted  signal  of  the  secondary  user  can  be 
expressed  as 


j- 1 

xs(t)  =  hWy/lhihit  -  iT)ej^t+^  (39) 

i= 0 

while  the  transmitted  signal  of  the  primary  user  is  given  by 

J- 1 

xp(t)  =  J2  &2(i)  V&2(*  -  iT)ej^t+^  (40) 

i=0 

where  Ek,  k  =  1,2  denotes  transmitted  energy  per  bit  and  <fikl  k  =  1,2  is  the 
carrier  phase  relative  to  the  transmitter’s  local  oscillator.  Information  bits  bk{i ), 
k  =  1,2  are  modulated  by  unique  digital  waveforms  dk{t),  k  =  1,2,  respectively 
given  by 


L  —  l 

dk(t)  =  T.  sk(l)gT{t  -  lTd),  k=  1,2  (41) 

2=0 

where  sjt(i)  £  {—1, 1}  is  the  l- th  waveform-bit  of  the  waveform  assigned  to  the 
fc-th  transmitter,  k  =  1,2,  gr{-)  is  a  pulse  shaping  square-root  raised  cosine 
(SRRC)  filter  of  duration  T^,  and  T  =  LT^  is  the  waveform  duration  (period). 

All  signals  are  considered  to  propagate  over  Rayleigh  multipath  fading  chan¬ 
nels  and  experience  complex  additive  white  Gaussian  noise  (AWGN)  at  the  re¬ 
ceiver.  Multipath  fading  is  modeled  by  a  linear  tappcd-delay  line  with  taps  that 
are  spaced  at  Td  intervals  and  are  weighted  by  independent  fading  coefficients. 
The  received  signal  in  the  presence  of  multipath  fading  is  bandlimited  to  1  /Td 
and  can  be  expressed  by 


N-l  N-l 

rc(t)  =  ^2  hi,nxs(t  -  nTd)  +  ^  h2,nxP{t  -  nTd)  +  n{t)  (42) 

n— 0  n— 0 

where  N  denotes  the  total  number  of  resolvable  paths  (assumed  to  be  the  same 
for  both  the  primary  and  secondary  transmitter),  k  =  1,  2,  n  =  0, . . . ,  N  — 
1,  are  independent  zero-mean  complex  Gaussian  random  variables  that  model 
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the  fading  phenomena  and  remain  constant  over  Tc.  n(t)  denotes  a  complex 
white  Gaussian  noise  process. 

We  note  that  low-cost  commodity  transceivers  introduce  carrier  frequency 
offsets  {A fk},  k  =  1,2,  between  any  transmitter-receiver  pair.  Therefore,  the 
received  baseband  signal  after  carrier  demodulation  is  given  by 

rW=EE^w 

2—0  k— 1 
N—l 

X  ^2  hk,ndk{t  —  iT  —  nTd)e~j2'KAtkt  +  n(t)  (43) 

n— 0 


where  hk,n  =  VEkhk^~^2n^cnTd~^k^  (the  resulting  carrier  phase  is  absorbed 
in  the  channel  coefficient).  After  pulse  matched  filtering  and  sampling  of  the 
received  signal  r(t)  at  rate  1  /Td  over  the  multipath  extended  symbol  period  of 
L  +  N  —  1  waveform-bits,  the  received  data 

r=[r(0)  r(Td)  ...  r((JL  +  Af-2)Td)]r  e  C(Ji+Ar_1)xl  (44) 

become  ready  for  filtering  and  detection  of  the  information  bits  of  the  secondary 
user  bi(i),  i  =  0, . . . ,  J  —  1.  With  respect  to  the  z-th  bit,  we  can  re-write  the 
received  data  vector  in  the  following  form 

y*  =  ir}(0+iL)Td-.Td:((L+N-2+iL)Td)  (45) 

which  is  equivalent  to 

Yi  =  b1(i)T1(i,  A/i)HiSi  +  p;  +ni;  i  =  0, . . . ,  J  —  1  (46) 

where  ri(0,A/i)  =  diag(l,  e-j27rAflTd, . . . ,  e-j(2^Afi(N+L-2)Td))  is  a  (jy  +  l  - 
l)x(N+L  —  l)  diagonal  matrix  consisting  of  the  frequency  offset  values  between 
the  secondary  transmitter  and  receiver,  and  Hi  £  q(l+n-i)xl  js  multipath 
fading  channel  matrix  given  as  follows 


N-l 

~h 1." 

n= 0 


OnxL 
I  LxL 

0(N-n-l)xL 


(47) 


We  assume  that  nj  ~  CAf(0,  ct^In+l-i)  represents  complex  zero-mean  white 
Gaussian  noise,  and  p,  accounts  for  the  unknown  primary  user  signal.  The 
frequency  offset  values  in  Ti(z,  A/i)  remain  constant  over  Tc  symbols. 


O.  Optimal  Waveform  Design 

The  cumulative  interference  for  the  secondary  user  signal  of  interest  in  (46)  is 
Pi  +  n.;.  If  A fi  and  the  channel  coefficients  hqo,  /ii  i, . . . ,  /iqiv-i  are  known 
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then  the  linear  filter  that  operates  on  y,  and  maximizes  the  SINR  at  its  output 
is  given  by 


WmaxSINR  = 


arg  maxw 


E{lwH(b1r1(i,  Af1)H1s1)|2} 
E{|wH(p  +  n)|2} 


=  +  (48) 


If  we  denote  the  cumulative  disturbance  autocorrelation  matrix  as  R/+at  =  Rp+ 
cr„I]v+i-i,  then  the  maximum  output  SINR  value  attained  at  the  secondary 
receiver  is  given  by 


SINRmax 


EllsTHfrH^.AfORj-^fbirRi.AfOHiS!)!2} 
E  {|STHHrH(i,  AfOR-^lp  +  n)|2} 


=  siTHfrf(i,A/i)R7H(ivri(j,A/i)H1si.  (49) 

Now,  we  consider  SINRmax  in  [76]  as  a  function  of  the  waveform  Si .  In  this 
project,  we  propose  to  design  s°pt  such  that  it  maximizes  the  SINR  at  the  output 
of  the  maximum  SINR  receiver  filter.  Then,  the  optimally  designed  waveform 
s°pt  is  fed  back  by  the  secondary  receiver  to  the  secondary  transmitter  enabling 
this  way  cognitive  channelization.  Thus,  s°pt  may  be  expressed  as 

s°pt  =  arg  maxSl  {siTHfrf(i,  AfijR^T^i,  Afi)HiSi}  .  (50) 

Let  us  now  define  R(A/i)  =  Afi)T{,J+Nri(i,  A/i)Hi.  Let  the 

eigenvectors  of  the  matrix  R(A fi)  are  denoted  by  qj ,  q2, . . . ,  qp  with  corre¬ 
sponding  eigenvalues  Ai  >  A2  >  •  •  •  >  Al,  then  s°pt  in  (50)  is  the  eigenvector 
that  corresponds  to  the  maximum  eigenvalue  Ai,  i.e. ,  s°pt  =  qi. 

We  note  that  waveform  optimization  under  the  constraint  of  maintaining 
the  performance  of  the  primary  user  above  a  presepecified  threshold  can  be 
carried  out  in  an  iterative/coupled  fashion  where  the  solution  provided  by  (50) 
is  used  to  further  optimize  the  power  allocated  to  the  secondary  user,  so  that  its 
interference  to  the  primary  user  is  maintained  below  a  prespecified  threshold. 
The  coupled  optimization  problem  will  be  addressed  in  a  continuation  of  this 
work,  and  is  not  considered  here. 


P.  Joint  Subspace-based  Channel  and  CFO  Estimation 

Expression  (50)  assumes  that  H1;  IbC  ■)  and  Ri+n  are  known.  In  practice,  we 
utilize  estimates  Hi,  Ti(-,  •)  and  R/+at  based  on  a  finite-size  of  data/observations. 
In  particular,  Rz+jv  is  sample-average  estimated  [82,  83]  over  J  silent  snapshots 
of  the  secondary  user,  i.e., 


R/+jv 


1 

J 


j 

y^(p,:  +  n,)(pj  +  nt)H. 

i=  1 


(51) 


47 


Approved  for  Public  Release;  Distribution  Unlimited 


The  rest  of  this  section  deals  with  the  problem  of  blind  estimation  of  the  un¬ 
known  channel  coefficients  and  frequency  offset  of  the  secondary  user  in  the 
presence  of  unknown  interference  and  AWG  noise.  The  phase  ambiguity  in¬ 
duced  by  blind  channel-estimation  methods  is  resolved  using  a  short  pilot  bit 
sequence  [80]. 

The  proposed  method  for  jointly  estimating  the  unknown  channel  coefficients 
hi>n  for  n  =  0, . . . ,  N  —  1  and  the  carrier  frequency  offset  A/1;  is  a  modified 
version  of  the  algorithm  presented  in  [80]  where  it  was  shown  that  the  rank  of  the 
noise  subspace  is  increased  by  two-sided  truncation  of  the  received  signal  vector 
y j.  For  our  setup,  the  later  implies  that  trimming  the  intersymbol  interference 
(ISI)  terms  for  the  secondary  user  from  both  sides,  reduces  the  signal  subspace 
rank  rs,  to  L  —  N+l  <  rs  <  2 K—  1  and  offers  the  maximum  possible  guaranteed 
minimum  rank  of  the  noise  subspace  of  (L  —  N  +  1)  —  (2 K  —  1). 

In  this  context,  let  the  (L  —  N  +  l)-length  truncated  version  of  the  received 
vector  y.(  with  respect  to  the  i-th  bit  be  as  follows 


yf  =  [Vi(N  -  1  )Td  Vi{NTd)  ...  yi((L-  1  )Td)] T  (52) 


which  is  equivalent  to 


(53) 


~  ~  ~  ~  T 

where  hi  =  [/ii.o  /ipi  . . .  /ii  jv— l]  is  the  channel  coefficients  vector  of  the 


secondary  user,  T*r(i,  A/i)  is  a  {L  —  N  +  1)  x  (L  —  N  +  1)  is  the  truncated 
version  of  the  diagonal  matrix  Ti(i,  A/i)  defined  in  (46),  and 


si[N-  1]  si [IV -2] 


si[0] 

*i[l] 


A  Sl[^]  s^N-1] 


(54) 


si[L  —  1]  s\[L  —  2] 


si[L-N] 


denotes  the  (L— N+l)  xN  matrix,  consisting  of  shifted  versions  of  the  waveform 
bits  of  the  secondary  user,  after  eliminating  N  —  1  elements  from  both  sides  of 
each  column. 


Let  Rtr  =  E{y[rylrH}  denote  the  autocorrelation  matrix  of  y[’r.  Let  Rtr 


QAQH  denote  the  eigendecomposition  of  the  input  autocorrelation  matrix, 
where  Q  is  a  matrix  with  columns  the  eigenvectors  of  Rtr  and  A  is  diago¬ 
nal  matrix  with  the  eigenvalues  of  Rtr-  Using  as  columns  the  eigenvectors  that 
correspond  to  the  (L  —  N  +  1)  —  (2 K  —  1)  smallest  eigenvalues,  we  form  the 
matrix  U„  £  £(l-~n+i)x((l-n+i)-(2K~i)) _  These  “bottom”  eigenvectors  of  Rtr 
belong  to  the  noise  subspace,  therefore  the  columns  of  U„  correspond  to  “noise 
eigenvectors” . 

We  note  that  when  Rtr  is  known,  the  nullspace  of  coincides  with  the 
signal  subspace  (ra  =  2 K  —  1).  However,  in  this  work,  Rtr  is  not  known, 
but  rather  sample-average  estimated  by  a  finite-size  record  of  truncated  input 
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data/observations,  i.e. , 


Rtr  =  7X>‘ryrH 

i= 1 


(55) 


where  J  denotes  the  size  of  the  data  record.  Thus,  U„  is  an  estimate  of  the 
noise  subspace,  based  on  the  eigendecomposition  of  the  estimated  input  auto¬ 
correlation  matrix  Rtr- 

The  estimates  of  hi  and  A/i  are  found  by  solving  the  following  constrained 
optimization  problem 

{Tf(i,Af1)As0h1)HtSn  =  0,  subject  to  ||hi||  =  1.  (56) 

The  solution  to  (56)  essentially  makes  the  secondary  user  signal  of  interest 
rf(i,  A/i)Aghi,  orthogonal  to  the  noise  subspace  [85].  Propositions  1  and  2 
below  are  modified  versions  of  Propositions  2  and  3  in  [80]  and  are  included 
here  for  easy  reference. 

Proposition  1  Expression  (56)  holds  true  for  values  of  A/i  chosen  to  make 
the  matrix 

Mx  (A/i)  =  Afrf^Af^V^Tf^Af^  (57) 

singular.  Given  the  choice  of  Afi,  we  select  hi  as  the  eigenvector  that  corre¬ 
sponds  to  the  zero-eigenvalue  o/Mi(A/i). 


The  following  proposition  provides  a  sufficient  condition  for  the  uniqueness 
of  the  solution  provided  in  Proposition  1. 

Proposition  2  The  parameters  Afi  and  hi  are  uniquely  determined  if 
(i)  rank(Ag)  =N  and 


(ii)  null(Ua )  D  range(rjr(i,  Af)Ag) 
_/  range(rf(i,Afi)As0h1) 
{0}, 


A  /  =  A/i 
A/  #  A/i. 


The  proposed  estimates  A/j  and  hi  are  then  used  to  form  wmaxsiNR,  and 
s^pt  as  described  in  Section  O.,  according  to  modified  versions  of  (48)  and  (50). 

We  note  that  the  blind  channel  estimate  hi  is  phase  ambiguous,  and  so 
is  wmaxsiNR-  Phase  ambiguity  is  resolved  by  applying  the  mean-square  where 
(f)  is  the  unknown  phase,  we  apply  the  mean-square  optimum  phase  recovery 
method  proposed  in  [80]  for  correcting  the  phase  of  the  overall  receiver  filter 
wmaxsiNR-  Specifically,  the  selection  criterion  of  f>  is  based  on  the  minimization 
of  the  mean-square  error  (MSE)  between  the  phase-corrected  filter  processed 
data  w^(axgINRyi  and  the  desired  information  bit  b\(i)  for  i  =  0, . . . ,  J  —  1, 

^  =  argmin0E{|w^axSINRyi-bi(i)|2},  <j>  €  [-w,ir).  (58) 
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The  optimum  phase  correction  according  to  the  optimization  criterion  (58)  is 
given  by  (see  [80]) 


=  angle  {wLxSINRE{yth}}  (59) 

where  E{yb][}  can  be  sample-average  estimated  by  j->  J2iLi  y ib*(i),  when  a  se¬ 
quence  of  P  pilot  information  bit  symbols  is  available.  Consequently,  we  recover 
the  phase  of  the  filter  estimate  wmaxgiNR,  then  we  correct  its  phase 


id> 

WmaxSINReJV, 


<j>  =  angle 


^niaxSINR 


1 

P 


^yib*(i) 


(60) 


t  L  i=l  J  ) 

and  finally,  we  apply  the  optimum  (in  terms  of  minimizing  probability  of  error) 
bit  detector  [86] 


Si(i)  =  sgn  |sR|w^axSINRe  J<Sq}}  ,  i  =  0, . . . ,  J  —  1.  (61) 


Linux  Host 


Fig.  12:  USRP2  block  diagram  with  attached  RFX2400. 


Q.  Experimental  Testbed  Setup 

In  this  section,  we  focus  on  the  software  design,  implementation  in  SDR  and  final 
deployment  of  a  low-cost  indoors  testbed  for  evaluating  the  concepts  of  cognitive 
channelization  presented  above.  For  this  reason  4  commercial  SDR  transceivers 
named  Universal  Software  Radio  Peripheral  (USRP2)  and  developed  by  Ettus 
[88]  are  used  to  provide  a  setup  of  both  a  primary  and  a  secondary  link. 

R.  USRP2  Receiver:  Low-cost,  Commercial  SDRs 

USRP2  is  a  low-cost,  simple  and  flexible  SDR  platform  that  consists  of  two  14- 
bit  analog-to-digital  converters  (ADCs),  capable  of  lOOMS/s,  two  16-bit  digital- 
to-analog  converters  (DACs),  capable  of  400MS/s,  a  Gigabit  Ethernet  (GigE) 
interface  and  a  Xilinx  Spartan  3A  FPGA  for  high  rate  signal  processing;  all 
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on-motherboard.  RF  front-end  functionality  is  provided  by  a  wide  range  of 
daughterboards  which  are  able  to  operate  from  DC  to  5.9GHz. 

In  our  experiments,  RFX2400  daughterboards  were  used.  These  boards  offer 
a  homodyne  (zero-IF  i.e. ,  no  intermediate  frequency)  receiver  architecture  and 
operate  at  2.25  —  2.9GHz  frequencies.  At  the  receiver  path  of  USRP2,  the 
FPGA  attached  on  the  motherboard  is  responsible  for  transferring  via  GigE  the 
incoming  from  RF  front-end  and  sampled  by  the  ADC,  digital  baseband  data 
to  a  host-PC  for  further  signal  processing.  In  the  case  that  USRP2  acts  as  a 
transmitter,  the  on-board  FPGA  is  waiting  for  PC  data  which  will  be  sent  to 
the  DAC,  get  converted  into  analog  samples  and  transmitted  via  the  attached 
daughterboard.  Received  baseband  signals  are  sent  to  the  PC  in  the  format  of 
16-bit  in-phase  and  16-bit  quadrature  data,  which  is  equivalent  to  4bytes  per 
complex  sample.  Therefore,  the  maximum  data  rate  for  the  baseband  IQ  data 
is  4B2™nIpie  —  30MS/s  (due  to  overhead  25MS/s).  However,  the  true  maximum 
sampling  rate  that  the  PC  is  capable  of  handling  depends  every  time  on  its 
hardware  specifications  and  data  processing  capabilities. 


Fig.  13:  Secondary  transmitter  block  diagram. 

Physical  layer  processing  takes  place  at  the  host  PC  by  using  GNU  Radio 
[89],  a  free  software  development  toolkit.  GNU  Radio  provides  the  signal  pro¬ 
cessing  runtime  and  processing  blocks  to  implement  software  radios  while  it 
controls  low-cost  external  RF  hardware  such  as  the  USRP2  through  the  Uni¬ 
versal  Hardware  Driver  (UHD).  UHD  provides  a  host  driver  and  an  application 
programming  interface  (API)  for  USRP2. 

GNU  Radio  offers  a  plethora  of  signal  processing  blocks  and  libraries  offer¬ 
ing  primitive  functionalities  i.e.,  filtering,  different  kinds  of  modulation  (BPSK, 
QPSK,  GMSK),  carrier  phase  synchronization  etc.  However,  due  to  the  unique 
and  novelty  features  characterizing  the  proposed  technique  of  cognitive  chan¬ 
nelization,  we  had  to  build  our  own  custom  signal  processing  blocks  and  wrap 
them  with  the  existing  interface.  Moreover,  the  timestamp  control  module  of¬ 
fered  along  with  GNU  Radio  is  exploited  in  our  scenario  for  enabling  accurate 
scheduling  between  the  secondary  receiver-transmitter  pair.  Specifically,  input 
and  output  streams  in  both  transceivers  are  tagged  with  timestamp  informa¬ 
tion,  allowing  us  in  that  way  to  control  the  silent  snapshots  of  the  secondary 
transmitter  during  which  the  sample-average  estimation  of  the  disturbance  au- 
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tocorrelation  matrix  R/+at  should  be  done.  Extreme  attention  was  also  given 
in  the  design  and  implementation  of  the  new  signal  processing  blocks  as  far 
as  their  number  of  input/output  streams  and  per  input  data  processing  is  con¬ 
cerned,  and  how  the  latter  is  correlated  with  the  throughput-delay  trade-off  that 
each  block  introduces,  and  overall  execution  time  of  the  block  in  the  runtime. 

S.  Transmitter  Design 

Both  transmitters  (secondary,  primary)  are  designed  and  implemented  in  order 
to  simulate  the  model  described  in  (39)  and  (40).  Fig.  13  depicts  the  block 
diagram  of  the  secondary  transmitter  implemented  in  GNU  Radio. 

First  a  random  byte  generator  source  feeds  the  transmitting  signal  process¬ 
ing  path  with  a  pre-specified  number  of  J / 8  bytes.  Then  each  byte  is  converted 
to  data  symbols  of  8  bits.  In  that  way,  we  come  up  with  a  sequence  of  J  bits 
£  {0,1},  where  each  of  them  is  mapped  to  either  1  or  —1  accordingly  to  the 
BPSK  constellation;  finally  we  get  £  { — 1,  -1-1}.  The  generated  bitstream 
is  then  given  as  an  input  to  our  own  custom  signal  processing  block  which  mod¬ 
ulates  each  bit  with  a  digital  waveform  of  bits  Sfc(Z),  l  =  0, . . .  L.  Moreover,  the 
incoming  stream  is  concatenated  with  a  stream  of  pilot  bits  who  play  the  role 
of  a  preamble  and  are  exploited  for  time  synchronization  and  phase  estimation- 
correction  at  the  receiver  side.  The  preamble  sequence  is  an  input  argument 
to  our  block,  while  the  waveform-bit  sequence  is  assumed  to  be  known  between 
every  transmitter-receiver  pair.  Waveforms  at  the  transmitter  are  dynamically 
changing  according  to  the  feedback  information  they  get  from  the  secondary  re¬ 
ceiver.  Finally,  pulse  shaping  is  applied  on  the  waveform-modulated  bitstream. 
For  this  reason  an  already  implemented  GNU  Radio  block  of  a  SRRC  filter  is 
used.  The  generated  pulse  train  at  the  output  of  the  last  block  provides  the 
in-phase  component  of  our  IQ  transmitter.  The  quadrature  component  is  given 
by  a  null  source. 

In  the  case  that  dk(l ),  k  =  1,  2,  l  =  1 ,L  is  fixed  and  not  changing  by 
sensing  the  feedback  information  of  the  secondary  receiver  the  design  proposed 
in  Fig.  13  fits  also  the  implementation  of  a  primary  transmitter. 


Fig.  14:  Secondary  receiver  block  diagram. 
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V.  Conclusions 


In  this  project,  we  proposed,  developed  and  studied  a  decentralized  algo¬ 
rithm  for  joint  dynamic  routing  and  code-division  channelization  in  cognitive 
radio  ad  hoc  networks.  We  considered  the  general  problem  of  maximizing  the 
network  throughput  through  joint  routing  and  spread-spectrum  channelization. 
We  proposed  an  algorithm  that  can  be  seen  as  a  distributed  localized  approxi¬ 
mation  of  the  throughput-maximizing  policy.  The  proposed  algorithm  requires 
solution  of  a  code-division  channelization  problem  as  the  search  for  the  sec¬ 
ondary  amplitude,  code  transmission  pair  that  maximizes  the  secondary  link 
output  SINR  subject  to  the  condition  that  all  primary  signal  output  SINR  val¬ 
ues  are  maintained  above  a  given  SINR-QoS  threshold  value.  The  formulated 
constrained  optimization  problem  is  non-convex  and  NP-hard  in  the  code  vector 
dimension.  Simulation  results  show  that  the  proposed  algorithm  considerably 
outperforms  baseline  solutions. 

We  also  presented  and  evaluated  through  experiments  RcUBe,  a  novel  radio 
framework  that  offers  abstractions  to  obtain  real-time  reconfigurability  and  self¬ 
optimization  capabilities  at  the  PHY,  MAC,  and  network  layers  of  the  protocol 
stack.  We  demonstrated  the  flexibility  of  RcUBe  by  implementing  a  cross-layer 
cognitive  radio  algorithm  designed  to  maximize  the  network  throughput  through 
jointly  optimal  control  of  routing,  dynamic  spectrum  allocation,  scheduling, 
transmit  power  control  and  interference  avoiding  waveform  selection. 
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List  of  Acronyms 


ADC 

AWGN 

BER 

BPSK 

CCC 

CSMA 

CDMA 

CDM 

CFO 

CTS 

CR 

CRN 

DAC 

DC 

DS -  CDMA 

DTS 

FDM 

FPGA 

ISI 

I  -  SMIRA 

LUT 

MAC 

MC -  CDMA 

MSE 

MSINR 

OFDM 

QCQP 

QPSK 

RTS 

ROCH 

RCA 

RcUBe 

SDR 

SINR 

SMI 

SRRC 

TDMA 

USRP 


analog  —  to  —  digital  converter 
additive  white  Gaussian  noise 
bit  —  error  —  rate 
binary  phase  —  shift  —  keying 
common  control  channel 
carrier  sense  multiple  access 
code  —  division  multiple  —  access 
code  —  division  multiplexing 
carrier  frequency  offset 
Clear  —  to  —  Send 
cognitive  radio 
cognitive  radio  network 
digital  —  to  —  analog  converter 
data  channel 

direct  —  sequence  code  —  division  multiple  —  access 

data  transmission  reservation 

frequency  —  division  multiplexing 

Field  —  Programmable  Gate  Arrays 

inter  —  symbol  —  interference 

interference  constraint  —  aware  stepwise  maximum  interference  removal  algorithm 
look  —  up  —  table 
medium  —  access  —  control 

multi  —  carrier  code  —  division  multiple  —  access 
mean  —  square  —  error 
maximum  —  SINR 

orthogonal  frequency  —  division  multiplexing 

quadratically  constrained  quadratic  program 

quadrature  phase  —  shift  —  keying 

Request  —  to  —  Send 

Routing  Code  —  division  Channelization 

random  code  assignment 

Real  —  time  Reconfigurable  Radio 

software  defined  radio 

signal  —  to  —  interference  —  plus  —  noise  ratio 
sample  —  matrix  —  inversion  method 
square  —  rootraisedcosine 
time  —  division  multiple  —  access 
Universal  Software  Radio  Peripheral 
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